Quantum kinetic approach for studying thermal transport in the presence of 
electron-electron interactions and disorder 
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A user friendly scheme based on the quantum kinetic equation is developed for studying thermal 
transport phenomena in the presence of interactions and disorder. We demonstrate that this 
scheme is suitable for both a systematic perturbative calculation as well as a general analysis. 
We believe that we present an adequate alternative to the Kubo formula, which for the thermal 
transport is rather cumbersome. 



Measurements of thermal and thermoelectric transport 
in electron systems of various kinds attracted consider- 
able attention in recent years.— Among such experiments 
one can find the measurement of the critical behavior of 
the thermopower near the metal insulator transition,^ the 
Nernst effect in high- Tc and conventional superconduc- 
tors^i^'^'^ thermal transport in quantum Hall systems^ 
as well as thermopower in quantum dots4 To benefit 
from these experimental efforts, theoretical studies are 
required. In this paper, we develop a new theoretical ap- 
paratus for analyzing thermal and thermoelectric trans- 
port in interacting electron systems in the presence of 
disorder. The strength of our scheme is in its generality 
that allow us to apply it for different kinds of interac- 
tions. We believe that this scheme can be adopted as an 
alternative to the Kubo formula. 

The validity of the Kubo formula for thermal and ther- 
moelectric conductivities was proved by Luttinger»2, A 
main ingredient in the Kubo formula is the quantum me- 
chanical expression for the current operators entering the 
correlation functions. The correlation functions describ- 
ing the thermoelectric or the thermal current involve the 
heat current operator which was derived by Luttinger^ 
for a system of interacting electrons. In the absence of 
the electron-electron interactions, the expression for the 
heat current operator is 



(1) 



+ -^^*™p(P'P')cj,(T)Cp'(T), 
P,P' 

where c^{t) (cp(t)) is the creation (annihilation) opera- 
tor of an electron in a state with energy Sp. Here /z is the 
chemical potential, Vimp{p,p') is the potential created 
by the disorder and r is the imaginary time. With the 
help of the equations of motion^" and after transforming 
to the Matsubara frequencies, the current operator can 
be written as: 



jft(q = 0,uj„) = 



dsp 2ien 



dp 



=p(e«)cp(e„ - ujn)- 
(2) 



When interactions between the electrons are included, 
the expression for j/i becomes more complicated. In gen- 
eral, even when the equations of motion are employed, 
the resulting expression for the current is not just the 
frequency multiplied by the velocity but rather contains 
additional termsj^ Unfortunately, the heat current op- 
erator of free electrons is frequently used for interacting 
systems even when there is no justification for it. We 
demonstrate in Appendix [B] that the simplified form of 
the Kubo formula fails to reproduce the phenomenolog- 
ical thermal conductivity of a Fermi- liquid system, i.e., 
the result docs not satisfy the Wiedemann-Franz law. 
The problems induced by the simplified Kubo formula do 
not necessarily imply that the use of the Kubo formula 
for thermal transport is generally wrong. The weak point 
is in replacing Luttinger's expression for the heat current 
by the simplified form. The problem with using the full 
expression for the heat current is that is that it is too 
complexjiS, 

The difficulties related to the Kubo formula lead us to 
turn to the quantum kinetic equation. Although deriva- 
tions of the transport coefficients using the kinetic equa- 
tion already exist, see for example Refs. ITslfl^flslllGlll?! . 
our method differs in few aspects. We developed a sim- 
ple scheme that, as we demonstrate in this paper, can 
be applied for various electron systems. The scheme is 
well suitable for both a general analysis and a systematic 
perturbation expansion. The novelty of our method is 
that we were able to derive both the kinetic equation and 
the currents directly from the action and the correspond- 
ing conservation laws. The systematic formulation of the 
quantum kinetic equation in the presence of a tempera- 
ture gradient has been achieved using Luttinger's idea of 
introducing a gravitational field.- As a result, we found 
the expressions for the currents from which all four com- 
ponents of the conductivity tensor can be extracted: 
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Moreover, we obtained that all the currents share a com- 
mon simple structure: 



de 
2^' 



dr'xe,hie) v(r,r',e)G(r',r,e) 
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Here we use the notation Xe(e) = — e for the electric cur- 
rent and Xhi^) — £ for the heat current. Both the vefocity 
and the quasiparticle Green's function in the above equa- 
tion are fully renormalized by the interaction. The cur- 
rents are related to the lesser component in the Keldysh 
space of the product of the renormalized velocity and 
Green's function. Equation 2] is one of the central re- 
sults of this paper; for the first time it is shown that 
the flow of energy occurs with the renormalized veloc- 
ity. This structure of the heat current guarantees that 
the Wiedemann-Franz law is satisfied in the framework 
of the Fermi-liquid theory. Finally, let us remark that al- 
though in this paper we consider electrons interacting via 
the Coulomb interaction, the generality of the method al- 
lows to account for different interactions with a minimal 
effort. In particular, this scheme was highly useful in the 
analysis of the Nernst effect in disordered films in the 
presence of superconducting fluctuations^^ 

The paper is organized as follows: In Sec. U we cal- 
culate the electric conductivity. We deliberately choose 
the well elaborated example of the electric conductivity 
in order to illustrate the main steps of our approach. We 
consider two scattering mechanisms; the Coulomb inter- 
action between the electrons and the elastic scattering by 
impurities. In Sec. [ill we repeat the scheme introduced 
in Sec. [Jin order to obtain the quantum kinetic equation 
for a system of electrons in the presence of a tempera- 
ture gradient. In Sec. IIIII we derive the expressions for 
the different currents induced by the temperature gradi- 
ent. For completeness, we present the expression for the 
heat current generated by an electric field. Sections [III 
and IIIII constitute the core of the paper. In the rest of 
the paper we demonstrate how to apply this scheme for 
various calculations and check that this method repro- 
duces some known results. In particular, in Sec. IIVI we 
show that the Wiedemann-Franz law is satisfied for a 
Fermi-liquid system. Then, in Sec.[V|we discuss the fate 
of the Wiedemann-Franz law when diffusive corrections 
arising due to the interplay of the electron-electron in- 
teraction and disorder are considered. In Sec. IVII we ex- 
amine Onsager's relations. The scheme developed in this 
paper allows us to calculate the two thermoelectric cur- 
rents separately. As an additional test, we demonstrate 
in Sec. lVII that our expressions for axx and otxx satisfy the 
Onsager relation, axx = Taxx- Appendix \K\ contains ad- 
ditional technical details of the derivation of the electric 
conductivity. In Appendix [Bj we present the calculation 
of the thermal conductivity for a Fermi-liquid system us- 
ing the simplified version of the Kubo formula. We show 
that the simplified formula leads to an erroneous result. 
In Appendix [Cl we concentrate on the contribution of the 
Coulomb drag to the different transport coefficients. We 
emphasize the difference between the role of the Coulomb 
drag in electric conductivity and that in thermal conduc- 
tivity. Finally, in Appendix [Ej we briefly describe the 
generalization of the scheme developed in this paper for 
superconducting fluctuations. 



I. THE QUANTUM KINETIC SCHEME FOR 
ANALYZING THE ELECTRIC CONDUCTIVITY 

In this paper we limit ourselves to the well known ex- 
ample of electrons interacting through the density chan- 
nel (e.g. the Coulomb interaction) in order to demon- 
strate the logic of our scheme for studying transport phe- 
nomena. For the electric conductivity we start with the 
following action: 



S= I dvdt \ i;\r,t)idti^{r,t) 



(VV't(r,t))(V^(r,t)) 



2m 



-^j''{r,t) [er •E + y,„p(r)]V'(r,t) 



1 / rfrVt(r,i)^t(r',t)f/(r-r')V(r',i)V'(r,t) 



(5) 



Here ip{r) and ip^'ir) are the Grasmannian fields describ- 
ing the quasiparticles, V'^(r) = X]p '-p^"*'^'"- electric 
field is assumed to be constant, and we choose to work 
in the gauge E — ~^ip{r). In general, the field oper- 
ators include a spin index that is summed over in the 
above action. Here and in the following we do not in- 
dicate (whenever it is possible) the spin indices because 
we do not consider any scattering mechanisms that flip 
the spins, or the Zeeman splitting. For convenience, we 
introduce the Hubbard-Stratonovich field (j){r) into the 
action: 

5 = / ,r,t { «.,r. 0.a*(r, «) - OKV^KM)) 

-^/.t (r, t) [er • E + F„„p(r)] V(r, t) - 0(r, t)V/ (r, i)iA(r, t) 



+ ^ I dv'<t>{v,t)U-\v~v')4>{v\t) 



(6) 



After this transformation the system is described 
by two propagators: the quasiparticle Green's func- 
tion, G{r,t;r' ,t'), and the propagator of the cj) fields, 
V{r,t;r' ,t'). [We use the term propagators when re- 
ferring to both these functions, while separately we 
name G(r, t; r', t') the quasiparticle Green's function, and 
V{r,t;r' ,t') the propagator of the interactions.] 

The general formulation of the quantum kinetic equa- 
tion for transport phenomena assumes that at time 
t = — oo the system is at equilibrium. Then, an ex- 
ternal field is adiabatically switched on generating cur- 
rents in the system. We follow the Keldysh- Schwinger ap- 
proach and use the matrix form of the propagators fi2i22iiSi 
We find it convenient to work in the basis of the re- 
tarded, advanced and Keldysh propagators. At equi- 
librium, the three components of the matrices are not 
independent, (e) = {1 - 2nF(e))(G^(e) - G^(e)) and 
V^{uj) = (l-h2np(a;))(y^(w)-y^(w)); herenF(e) is the 
Fermi distribution function and rtp(a;) is the Bose distri- 
bution function. The equilibrium propagators give the 
properties of the system at t = —oo. In order to describe 
the system in the non-equilibrium state caused by the ex- 
ternal field, one should derive the system of equations for 
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FIG. 1: (a) Illustration of Eq.[8a]for the full Green's function G. 
(b) The Dyson equation for Gi„t (see Eg. I8bp . Note that Gint 
includes scattering by impurities only through EjC, V}, which 
is a functional of the propagators G and V. The bare Green's 
function, denoted by Gb, is free from the interactions and the 
scattering by impurities. 



G'i„t(r, r', t') includes partially the scattering by im- 
purities. We will rewrite the Green's function and the 
kinetic equation in terms of the center of mass coordi- 
nates, R = (r + r')/2, T = (t + t')/2, and the relative 
coordinates, p = r — r' , t = t — t' . There are two sources 
for the dependence of the Green's function on the center 
of mass coordinate, the electric field and the scattering 
by the impurities. To separate these two dependencies, 
we use the notation G(R, T; p, r; imp), where the explicit 
dependence on R is caused by the electric field, while the 
one arising due to the impurity potential is incorporated 
into imp. Next, we Fourier transform the Green's func- 
tion with respect to the relative time coordinate. Simul- 
taneously, we wish to write the quantum kinetic equa- 
tion in the gauge invariant form. Owing to the fact that 
the gauge invariant time derivative includes the scalar 
potential, one can modify the Fourier transform in the 
following wayj ^^'^^ 



the propagators, i.e., the quantum kinetic equation. The 
starting point is the Dyson equations corresponding to 
the action which includes the external field. The Dyson 
equation for the electric field dependent Green's function 
of the quasiparticles is: 



G{v,t-v\t') (7) 
5{v - Y')5{t -t') + l dhdr,j:{r, t; r„t,)Gir„U;r' ,t'). 



The chemical potential appears because we work in the 
grand canonical ensemble. Note that the specific time de- 
pendence of the retarded, advanced and Keldysh propa- 
gators is encoded in the definition of the different compo- 
nents of the matrix propagators. Correspondingly, when 
the matrix propagator is Fourier transformed, each of its 
components acquires the proper analytic structure in the 
complex frequency plane. 

We choose to postpone the averaging over impuri- 
ties until the last stage of the derivation of the current. 
Therefore, the Green's function of the quasiparticles con- 
tains open impurity lines as illustrated in the two coupled 
equations presented in Fig. [TJ 

G(r,t;r',i') = G,„t(r,i;r',i') (8a) 

+ / dri(itiGi„f(r,t;ri,ti)T/i,„p(ri)G(ri,ii;r',t'); 



G(R, T; p, e; imp) = J dre'^'+'^^^^GiR, T; p, t; imp). 

(9) 



The underscore is used to mark gauge invariant quanti- 
ties. 

After linearizing the quantum kinetic equation for G 
with respect to the electric field, one obtains: 



>^ Q.{p,e;imp) — 6{p) + / driS (p — ri, e; imp) (10) 



eE 



d . .d 



G (ri, e; imp) . 



Gint{r,t;r' ,t') = Gh(r, i; r', t') +jdVidtidv2dUGb{r,t]T^,ti) 
X t{Y„t,-v„t,)G,nt{v„U-y,t'). (8b) 

Here, Gmi(r, i; r', t') is the Green's function of inter- 
acting electrons, while Gb(r, r', t') is free from both 
interactions and scattering by impurities. Note, that 



A detailed derivation is presented in Appendix lAl Here, 
V = |Vfl -I- Vp. Clearly, for a constant electric field 
there is no explicit dependence of the gauge invariant 
Green's function on the center of mass coordinate, and 
the entire dependence of G(p, e; imp) on R is incorpo- 
rated into imp. To find G, we write the Green's function 
as a sum of two terms; the equilibrium Green's function 
geq{p,e;imp) and the E-dependent part G-EipT^^imp). 
The equation for the equilibrium Green's function is: 



e+^+M-^imp geqip,e;imp) ^ S{p) (11) 
- J dr^&eqip - Ti, e; imp)geq{ri,e; imp). 



Notice that the Laplacian in the above equation includes 
the derivatives with respect to the center of mass and 
relative coordinates. [The dependence on the center of 
mass coordinate is caused by Vimp.] At equilibrium, the 
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components of the matrix Green's function are: 



gf/(p,e;zmp) (12) 
= [e + VV2m + n~ V,^p - crf/(p, e; imp)] 

9eq{p, e; imp) 

= (1 - 2nF{e))[gfq{p,e\imp) - g^q{p,e;imp)]. 

Here gi^^'^' is analytic on the upper (lower) half of the 
complex frequency plane e. 

The detailed calculation of Ge presented in Ap- 
pendix]^ yields : 



Ge(p, e) = geq (e) (e) 5eg (e) 

ieE (e) , .^ » /^- t \ ^9eq (e) 
Veg(e)5e9 (e) - 5eg (e) Ve,(e) 



(13) 



5e 



From now on, whenever the dependence on the coordi- 
nates in the product of matrices is not specified, it should 
be understood as a convolution in real space. The matrix 
is the velocity of the quasiparticles renormalized by 
the interactions at equilibrium: 



V - V' 

Veg(r,r',e) = -i lim — i(r - r')CTe,(r, r', e) 

r'-*r 2m 



(14) 



The equation for Ge contains the electric field dependent 
self-energy, Se, which by itself is a function of the E- 
dependent propagators. Thus, in order to find Ge, one 
has to determine the structure of the self-energy. Once 
the form of the self-energy is fixed, one should take into 
consideration that each of the propagators in Ee may 
depend on the electric field. 

To find the response to the applied electric field, we 
need also to derive the quantum kinetic equation for the 
propagator of the interaction. The Dyson equation for V 
is: 



dr,dt,U-'^{r - r,)6{t - t,)V{r„t,;r' , t') = S{r - r') 
dr,dtA{r, t; r,,t,)V{r„t,;r' , t'), (15) 



where 11 is the self-energy (polarization operator) for the 
interaction field (j){r,t). Unlike the kinetic equation for 
the quasiparticle Green's function, the electric field does 
not appear explicitly in the equation for V. Indeed, since 
the interaction field is neutral, the electric field can 
only enter through the self-energy term 11 that contains 
the Green's functions of the charged quasiparticles G. 
For the same reason, the Fourier transform of the rel- 
ative time in V{r,t;r' ,t') should be performed without 
the gauge factor. 

The calculation of the electric conductivity requires 
the expression for the electric current in terms of the 
propagators depending on the electric field. We derive 
the electric current through the continuity equation. We 



start from the density which can be related to the lesser 
Green's function in the following way: 



n{r,t)^2 lim {ip'' {r' ,t')'il;{r,t)) (16) 

r'— *r 

= -2i lim G<(r,i;r',i')- 



We use the notation t' i+ to indicate that the limit 
should be taken in such a way that t is on the upper 
branch of the Keldysh contour, while t' is on the lower 
branchi^ The summation over the spin projection results 
in a factor of 2. The lesser component of the Green's 
function can be written in terms of the retarded, ad- 
vanced and Keldysh Green's functions through the rela- 
tion G< = (G^ -G^ + G^)/2. Here {A) denotes the 
quantum mechanical averaging with the action given in 
Eq. [6] Therefore, the Green's function is fully dressed by 
the interactions and depends on the impurity potential. 
In addition, G is a function of the electric field. Since 
we find the current by extracting it from the continuity 
equation, we assume that the electric field has some spa- 
tial modulations that will be set to zero at the end of the 
procedure. 

The continuity equation for the charge density, 
— eri(r, t) + Vje(r, t) = 0, can be be written as: 



Vje(r,t)-2e hm 1 + A ) (v/(r', OV'(r, i)) 



(17) 



Under the average, the equations of motion for the fields 
ip and ip'^ allow us to rewrite the continuity equation as 
follows: 



Vje — 2ie lim 

r' — 'V 

t'-,t+ 



2m 



e(r - r')E - V,mp{v) 



+V^^™p(r')-'/'(r,t)+0(r',t') 



^t(r',i')V'(r,t) 



(18) 



We may express the RHS of the above equation in terms 
of the Green's functions and self-energies: 



Vje = 2e lim 



y2 _ y/2 

2m 



G(r,i;r',t) 



(19) 



dv^dtj:{v, t; Ti, ii)G(ri, ii; r', t) 



dr,dt,G{r, t; r,,t,)t{r,,t,;r' , t) 



We use the notation [...] ^ to indicate that the expression 
inside the square brackets is a matrix and the current 
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corresponds to the lesser component of this matrix. The 
exphcit dependence on E in Eq. [TH] dropped out as a 
result of taking the limit r' r. However, the Green's 
functions as well as the self-energies in the expression for 
the electric current depend on the electric field. 

To resolve the expression for the electric current, we 
have to represent the RHS of Eq.[Tn]as a gradient of some 
function. With this in mind, we write the convolutions 
of the Green's functions and self-energies in terms of the 
center of mass and relative coordinates: 



GE - EG 



lim / dvidt. 



G ( R + ^^—2^^; ^ - Ti, -ti]imp 



xG IR+ ''' ^^^'^ ;^+r,,t,;imp 



(20) 



To extract the current out of the continuity equation we 
expand [EG — GE] ^ with respect to the deviation from 
the center of mass coordinate. [This step resembles the 
gradient expansion discussed, for example, in Refs. [20ll2lj | 
As a consequence of the structure of Eq. [20l one may im- 
mediately see that all even terms in the expansion vanish. 
Generally speaking, the expansion includes both the ex- 
plicit dependence on R arising due to the external field 
and the one entering through the impurity potential. Re- 



call that ultimately we are interested in the current aver- 
aged over space. Owing to the fact that any correlation 
function between the impurity centers depends only on 
their relative distance, the derivatives of Vimp with re- 
spect to the center of mass coordinate vanish upon aver- 
aging. Therefore, in the regime of linear response, it is 
enough to keep only the first non-vanishing term in the 
expansion: 



GE - EG 



(21) 



XyE (R; ri,ti;imp) 



h.c. 



After performing the gauge invariant Fourier transform 
on the current, one come to the very compact expression: 



le I —dv' 

2-K 



v(r,r',e)G(r',r,e) 



-I- h.c. 



(22) 



We used the fact that G and v in the expression for the 
current have the same center of mass coordinate, while 
their relative time coordinates differ by a sign. There- 
fore, when we Fourier transform both the velocity and the 
Green's function according to Eq. [HI the two gauge fac- 
tors cancel each other. As a result, je becomes a simple 
convolution of two gauge invariant quantities. The ex- 
pression for the current contains the renormalized veloc- 
ity as defined in Eq. [H] with the only difference that the 
self-energy may depend on the electric field, v = Veg+VE, 
where VE(r, r', e) = —i{v — r')EE(i", i"', e)- 



The final step in the derivation of the electric conductivity is to insert the expressions for Ge and ve into Eq. [22l 
and complete linearizing it with respect to the electric field (for more detailed see the end of Appendix |Aj . Using the 
known relations between the components of the Green's function at equilibrium, one gets: 



27r dt 



-g^,{e)vfie)g^,{e)vf{e)]-e'E, / -np{e) 



de 



«f (.)^.f (e)4(e) - vf{e)g^^{e)vfie)^ 



de 



de 



/ ^«f(e)4(e) [S<(e)(l - UFie)) + E>(6)np(6)] {g^^{e) - g^^{e)) + c.c. 



The second term in the above expression is zero for the 
longitudinal current, but it becomes important when the 
transverse current in the presence of a magnetic field is 
considered. 

The obtained result for the electric current is expressed 
in terms of the renormalized Green's functions and veloc- 
ities, i.e., it holds to all order with respect to the electron- 
electron interaction. In practice, to get a quantitative 
result, one has to specify the form of the self-energy and 
to average over the disorder. In the case of the electric 



conductivity, the terms in Eq. [23] are equivalent to those 
given by the Kubo formula after performing in the lat- 
ter the analytic continuation to the real frequency. The 
derivatives with respect to the frequency are the same 
derivatives that one gets after expanding the Kubo for- 
mula with respect to the external frequency. While the 
derivatives with respect to the frequency in the first two 
terms of Ea. l23l appear explicitly, in the last term they re- 
veal themselves only after linearizing the self-energy with 
respect to E. If interested, one can generate the pertur- 
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bative expansion order by order in a systematic fashion, 
and give a diagrammatic interpretation for each of the 
terms. This technique automatically determines the an- 
alytic structure, the way the distribution functions enter 
as well as the numerical factors of all the diagrams. [We 
checked that for the electric conductivity the set of di- 
agrams obtained in the quantum kinetic approach coin- 
cides with the one given by the Kubo formula. We re- 
produced all the diagrams contributing to the Altshuler- 
Aronov corrections to the conductivity)^!^ In order to 
get these known corrections, one should properly perform 
the averaging over the disorder in Eq. [531] 

To summarize, we demonstrated the main steps in the 
derivation of currents in response to an external field us- 
ing the electric conductivity as an instructive example. 
In particular, we showed how to find from the continu- 
ity equation the current in terms of the renormalized 
(dressed) quantities. In the subsequent sections we shall 
follow the same scheme in order to find the heat and elec- 
tric currents as a response to a temperature gradient. In 
addition, in Appendix |E] we describe how this scheme can 
be applied for an interaction in the Cooper channel. 



II. DERIVATION OF THE KINETIC 
EQUATION IN THE PRESENCE OF A 
TEMPERATURE GRADIENT 

Similar to the calculation of the electric conductivity 
described in the previous section, the derivation of the 
electric and heat currents as a response to a tempera- 
ture gradient consists of two steps. One has to derive 
the kinetic equation and to find the expressions for the 
currents from the continuity equations. 

We start with the quantum kinetic equation for the 
matrix Green's function. Since the temperature gradi- 
ent is not a mechanical force, one cannot obtain the re- 
sponse to VT just following the route we elaborated for 
the electric field. To overcome this obstacle, we (follow- 
ing Luttinger— ) introduce an auxiliary gravitational field 
that enters the action. We will show how to establish 
a direct connection between the response to the gravita- 
tional field and the response to the temperature gradient. 

In the Keldysh-Schwinger approachfi2iS£i2i the Green's 
functions arc defined using the time ordering operator Tc 
along the Kcldysh contour C (see Fig. [2]): 

G(r,t;r',i') (24) 
= ~i (Tc [e-'^o drmT)-^N)^(^^^ t')})- 

The integration in the exponent contains two parts. The 
first part of the contour is parallel to the imaginary axis 
starting at t = —oo + i/3 and ending at t — ~oo + iS. The 
second, which is parallel to the real time axis, gives the 
evolution of the system in time. The Green's function 
along the second part of the contour is described by the 
kinetic equation. In the derivation of transport proper- 
ties we usually assume that the driving force is switched 



Imt 



t=-oo+ii: 



t-Hi5 



t-iS 



FIG. 2: The Keldysh contour. The first, vertical, part of the 
contour is parallel to the imaginary axis starting a.t t = —oo + ifi 
and ending at t — —oo + iS. The second, horizontal, part of the 
contour is parallel to the real time axis. Here /3 is the inverse 
temperature. 



on adiabatically, starting at t = —oo + iS. Thus, the in- 
tegration of the Hamiltonian in the exponent along the 
first part of the contour yields the thermal distribution 
of a system in the unperturbed state. As a result, at 
t — —oo + iS when the external perturbation starts to 
act, the system is at thermal equilibrium. 

In contrast to the electric field described in the preced- 
ing section, a space dependent temperature, T(r), infiu- 
ences the first part of the integration. In principle, one 
may try to generalize the integration along the first part 
of the Keldysh contour in the following way: 

~i dr dt {h{r) - ^n{r)) \ , (25) 

J J-oo+i/3(r) J 

where h(r) is the Hamiltonian density. Clearly, the den- 
sity matrix obtained as a result of this integration corre- 
sponds to a non-uniform state. On the other hand, the 
integration along the second part of the Keldysh contour 
(parallel to the real axis) is independent of the temper- 
ature and, therefore, the quantum kinetic equation does 
not include any external perturbation. In other words, 
we face the problem of finding highly non-trivial initial 
state before we even start to study its time evolution. 
To avoid this complicated task, we shall reformulate the 
problem in such a way that the initial state of the system 
is uniform in space. 

A similar problem of treating a spatial varying initial 
state appears when the response to a gradient of the den- 
sity is studied. According to Einstein's construction, in 
the case of a density gradient the stationary state can be 
obtained by adding a scalar potential at time t = —oo 
in such a way that the electro-chemical potential is kept 
constant, C(r) = /i(r) — e(p(r) — const, and the initial 
state is uniform. Then, the response to the gradient of 
the chemical potential (density) can be derived by adi- 
abatically switching off the scalar potential, ip{r). [This 
is one way to interpret the Einstein relation.] In a sim- 
ilar fashion, Luttinger^ introduced a gravitational field, 
7(r), as a counterpart of the non- uniform temperature. 
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This auxiliary field allows deriving the equivalent of the 
Einstein relation for studying thermal transport. 

We wish to introduce Luttinger's gravitational field 
7(r) into the quantum kinetic approach. Since the grav- 
itational field can be considered as a spatial dependent 
measure of the time coordinate, we shall study the fol- 
lowing action: 



(26) 



(V7At(r,t))(VV(r,t)) 



2m 

-0(r, t) (r, t) V(r, t) + ^ J dv'<t>{v, t)U-\v - r')</)(r', t) 

Let us first fix 7(r) in such a way that under the combined 
effect of the spatially dependent temperature, T(r) = T+ 
ST{r), and the gravitational field the system remains at 
equilibrium with a uniform effective temperature jqT. In 
other words, the distribution function that describes the 
state of the system does not evolve in time, and it is equal 
to nF{e) = [e<=/(To^) + ■ One may notice that the 
gravitational field enters Ea.[26lonlv through the product 
t"f{r). Correspondingly, after Fourier transforming the 
time coordinate, this field is coupled to the frequency 
as e/j{r). Therefore, under the combined effect of the 
temperature and the gravitational field, the equilibrium 
distribution function becomes 



exp 



7(r)T(r) 



(27) 



where, the gravitational field has to be chosen so that 
the product 7(r) • T(r) — 70T is constant in space. [We 
keep 7o unspecified until the end of the derivation of the 
thermal and thermoelectric currents when it is set to be 
1. We shall see that our choice to leave 70 simplifies the 
analysis.] 

Next, in order to find the response of the system to 
the temperature gradient, we adiabatically switch off the 
gravitational field starting at time t = — 00 + iS (or equiv- 
alently switch on the same field with an opposite sign). 
Now, the situation is similar to the one we encountered in 
studying electric conductivity. While for the first part of 
the integration, parallel to the imaginary axis, the com- 
bined effect of the temperature and the gravitational field 
maintains the system at equilibrium, in the second part 
of the contour integration the change in the gravitational 
field perturbs the system. 

In the following we shall study the kinetic equation in 
the presence of the gravitational field as described by the 
action in Eq. 1261 Since we are specifically interested in 
switching off the field 7(r) — 7oT/r(r), in the end of 
the derivation we will set <57(r) = ST{r)/T. [The last 
relation holds for the linear response. Pay attention to 
the sign in this relation. It is a consequence of the fact 
that switching off the field 6^{r) is equivalent to switch- 
ing on the field —Sj{r).] According to Eq.l^H the Dyson 



equation for the Green's function in the presence of the 
gravitational field is: 



d V(7(r)V) 

^^ + ^^-7(r)(Wr)-.) 



= 6{r - r')6{t - t') 
+ 7(r) / dtidriE(r,i;ri,i07(i-i)G(ri,ti;r',t') 



G(r,i;r',0 
(28) 



The collision integral contains 7(ri) because each inte- 
gration over time is accompanied by the gravitational 
field. At this stage, if to write the gravitational field 
as 7(r) = 7o -f rV7, it is possible to expand straight- 
forwardly the kinetic equation up to the linear order in 
V7. However, a radical simplification may be achieved 
by applying the following transformation to the kinetic 
equation: 

y(r, t; r', t') = 7-'/'(r)£(r, t; r', t')r^/\v'), (29) 

where Y can be the Green's function or the self-energy. 
When terms of the order (V7)^ and V^7 are neglected, 
the equation for the Green's function G(r, i;r',t') ac- 
quires the form: 



d 



"7(r)9t ' 2^'^"''''^''^^ 
X Gir,t;r' ,t') ^ S{r - r')S{t ~ t') 



'1/2, 



r') 



(30) 



+ 7i/2(r) / dt,dv,t{r,t;r„t,)G{r„t,;r',t')^-^^^{r'). 



Multiplying the above equation by 7~^/^(r) and 7^/^(r'), 
we eliminate the gravitational field from all terms beside 
the time derivative. 

After Fourier transforming the relative time in Eq. \30\ 
the kinetic equation in the regime of linear response be- 
comes: 



1 



RV7 pV7\ e 



70 



70 / 70 



2m 



-F™p(R + p/2) + A* > G(R; p, e) = S{p) 



R 



(31) 



J dr,t(^+^;p-r,,e)a 



Here we separated the dependence on the center of mass 
and the relative coordinates. Once again, we used V = 
+ Vp. Since the gravitational field is independent 
of time, and owing to the fact that we are interested 
in the stationary state, we may omit the center of mass 
time. It is worth paying attention to the difference in the 
way the electric and gravitational fields enter the kinetic 
equation. The electric field in the gauge invariant kinetic 
equation appears only with the relative coordinate p. As 
a result, the gauge invariant Green's function depends on 
the center of mass coordinate only due to the scattering 
by the impurities. [Recall that we postpone the averaging 
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over the disorder till the end of the calculation.] The 
essential feature of the quantum kinetic equation in the 
presence of a gravitational field is that V7 multiplies not 
only p but also R. Consequently, the gravitational field 
induces an additional dependence of the Green's function 
on the center of mass coordinate. 

We shall write the Green's function G as a sum of the 
equilibrium Green's function g^q and the V7-dependent 
Green's function. The kinetic equation describing g^q can 
be obtained by setting V7 — 0: 



— + 7; V;™p(R + p/2) +^ 

70 2m 



geq{p,e]imp) (32) 
= 5(p) + j dv^deq (p - ri,e]imp)geq {r^,e]imp) . 



The above equation practicably coincides with Eq. 1111 
the only difference is that here the frequency appears as 
e/70. As we have already done in the preceding section, 
we separate the dependence on R due to the driving force 
from the one caused by the impurity potential. In our no- 
tation, the latter has been incorporated into imp. Before 
averaging, the dependence of the Green's function on the 



disorder is described by the set of equations which was 
already given in Eq.[51 After averaging over the disorder, 
the equilibrium Green's function becomes: 



{9^q\P.^))., 



(33) 



2m 



6))_^ = (1 - 2nF (6/70)) [4(p, 6) - 4(p, 6)] 



where r is the mean free time. These Green's functions 
can be interpreted as describing the quasiparticles in the 
equilibrium state with temperature T, and effective fre- 
quency e/70. 

As we already discussed, the explicit dependence on 
R in G(R; p, e; , imp) is induced by the gradient of the 
gravitational field. Therefore, in the process of linearizing 
the kinetic equation with respect to V7, we expand G 

and E in the collision integral with respect to this explicit 
dependence on R. In other words, we may rewrite the 
last term in Eq. [5l1as: 



drig (R + y ;P - ri,e;imp^ ^ ( ^ e; ^'^P 



Ti 9S(R; p - Ti, e; imp) 



9R 



G(R; Ti , e; imp) — j drj S(R; p — rj , e; imp) 



(iriS(R; p — Tj, e; imp)G(R; ri, e; imp) 
p - Ti dQ_{Ii; Ti , e; imp) 



9R 



(34) 



We will see that the last two terms are indeed propor- 
tional to V7. 

We separate the part of G depending on V7 into two 
pieces. The equation for the first one, Gioc-eq, is: 
72 



IS 



70 2m 

RV7 

—<^9eq{p, e; imp) 

7o 



Gzoc-eg(R; P, e; imp) 



(35) 



= / dY^aeq{p-r^,e;imp)Gioc-eq{^]Y^,c,imp) 
dri t^ioc-eq (R; p - r 1 , e; imp) g^q (ri , e; im,p) . 



Following the steps presented in Appendix |^ for the 
derivation of the E-dependent Green's function, we may 
rewrite the expression for Gioc-eq (R; P, e; imp) as: 



Gloc-eq{() = (RV7)ge,(e)— ge?(e) 

To 

+ (R-V7)ge?(e)Sioc-eg(e)5eg(e) 



(36) 



Once again, one should understand the product as a con- 
volution of the coordinates. The solution of this equation 



Gioc-eq (R; P, e; imp) = (RV7) 



]ioc-eg(R; P, e; imp) = (RV7) 



dgeqjp, e; imp) ^ 

dd-eq{p,f^;imp) 
970 



(37a) 



(37b) 



We see that the local equilibrium Green's function, 
Gioc-eq, is a straightforward extension of the equilib- 
rium Green's function for a non-uniform gravitational 
field/temperature. This part of the Green's function 
describes the readjustment of quasiparticles to the non- 
uniform gravitational field/temperature when the system 
is trying to maintain a local equilibrium. This response 
of the electrons to the non-uniform perturbation tempts 
to induce a spatial modulation of the density. Since for 
charged particles it is impossible to have a large scale 
charge modulation, the gradient of the gravitational field 
transfers into a gradient of the electro-chemical potential. 
Therefore, je = (t(E -|- V/x/e) — ctE* where the effective 
field E* is the one measured in experiments. In other 
words, although this contribution to the current is initi- 
ated by the temperature gradient, it reveals itself through 
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the electric conductivity. We wish to remark that for a 
constant electric field an equivalent for the local equi- 
librium Green's function does not appear. [The role of 
the local-equilibrium Green's function is most peculiar 
when the response to V7 is considered in the presence 
of a magnetic field. Under these conditions, Gioc-eq is 
responsible for the non-vanishing contribution to je from 



the magnetization current 1^] 

The equation for the second term of the V7- 
dependent part of the Green's function, denoted by 
G-^'y(R; p,e;ijnp), includes all the terms in Eg. 1311 that 
did not enter the equations for geq and Gioc-eq- 



p-r,,e; ^mp) Gv, (R; r„ e; ^mp) - ^eg^^ e; ^mp) + ' ^^^.„._.,(R; p, ^mp) 



2m 



dRdp 



(38) 



driSv7 (R-; P — Ti, e; imp) (ri, e; + / dr^ 



Ti dJ^ioc-eq (R; p - Tj, e; imp) 



9R 



geq (ri,e;imp) 



p-r^ dGioc-eq (R; fi, e; imp) 

2 on 



In the above equation, the derivatives with respect to the 
center of mass coordinate act only on the explicit depen- 
dence of G/oc-eg(R, P, e;im-?3) and E;oc-eg(R, P, e; ^"^^') 
on R (i.e., through the spatial dependent gravitational 
field). Note that the derivatives with respect to the cen- 
ter of mass coordinate which act on Vimp in the local equi- 
librium Green's function was already included in Eq. 1351 
Once the explicit expressions for Gioc-eq and Ejoc-eq 
are inserted into Eq. I38[ the kinetic equation for G V7 be- 
comes similar to the one for Ge analyzed in Appendix. |A] 
(see Eq. IA2p . After similar manipulations, the kinetic 
equation acquires a simple form resembling Eq. 1131 



Gv7(p, e; , imp) = g^q (e) t^^ (e) g^q (e) 



(39) 



iV7 



dgeq (e) 



Veq{e)geq (e) - 5eg (e) Veg(e)- 



Jeq 



The expression for the renormalized velocity was already 
defined in Eq. [TH Let us emphasize that despite of the 
similarity in the structure, the equations for Ge and Gv7 
are not identical. The derivative with respect to the fre- 
quency in Eq. [13] is replaced by a derivative with respect 
to 7o in Eq. [351 Since the derivative with respect to 70 
acts on the equilibrium Green's functions for which 70 
accompanies the frequency as e/70, this derivative can 
be replaced by — (e/7g)9/9e. We see that according to 
the quantum kinetic equation for Gv7, the derivative 
with respect to the frequency is multiplied by the same 
frequency. As we have pointed out in Sec. [H the deriva- 
tive d I de corresponds to the expansion with respect to 
the external frequency in the Kubo formula. In the sim- 
plified version of the Kubo formula given by Eq. [21 the 
frequency multiplies the unrenormalized velocity vertex, 
which is not connected to the expansion with respect to 
the external frequency. As a result, in the presence of 
interactions, the frequency in the derivative and the one 
multiplying the velocity vertex are not necessarily the 



same. The difference between the expressions obtained 
by using the two methods may seem minor. However, it 
will become clear that this subtle point actually leads to 
different expressions for the thermal conductivity even in 
the Fermi-liquid theory. 

To complete the derivation of the response of the sys- 
tem to V7, we have to find the dependence of the prop- 
agator of the interaction V on the gravitational field. 
Let us return to the action presented in Eq. 1261 In 
general, for a non-local in space interaction, one should 
be concerned about which of the coordinates, r or r', 
should be prescribed to 7. This question is relevant 
for the last term in the action which describes the bare 
interaction via the field (p. In this paper we consider 
electrons interacting only through the Coulomb interac- 
tion. If one recalls that for the Coulomb interaction, 
/ dr'0(r,t)f/-i(r - r')0(r',i) = e2(V0(r, t)) VStt, it is 
clear that the problem of attributing the coordinate to 7 
does not exist. 

The Dyson equation for the propagator of the interac- 
tion in the presence of a gravitational field is: 



e 

8^ 



V (7(r) (r, t; r', t')) = 5{y - v')5{t - t') (40) 
7(r) / dv,dtJl{v,t-v,,t,)-i{v,)V{v,,t,-v',t'). 



In the derivation of the above equation we used the spe- 
cific expression for the Coulomb interaction. An impor- 
tant feature of the kinetic equation for the propagator 
of the instantaneous interaction is that it does not in- 
clude any time derivatives. As a consequence of this fact, 
the explicit dependence on the gravitational field may be 
eliminated from the equation in the linear response by 
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transforming to the propagator V_: 

-J dr,dtMr,t;r„t,)-g{r,,t,;r\t'), (41) 

Here we employ the same transformation as in Eq. 1291 
As one can see, the entire dependence of the propagator 
V_ on the gravitational field is through the quasiparti- 

cle Green's functions that enter the self-energy ij. Note 
the similarity between the kinetic equation for the prop- 
agator of the interaction in the density channel in the 
presence of a gravitational field and Eq. [13 

Let us separate the solution of Eq. 21] into the equilib- 
rium and V7-dependent propagators, y_ — Veq+Vioc-eq + 
Vv7- The propagator at equilibrium satisfies the equa- 
tion: 



Veqi-R;p,Lo)^ U-\p) +n,qCR; p,Lo) 



(42) 



The entire dependence of Veq (R; p, uj) on the frequency 
is due to the quasiparticle Green's functions in 11. Hence, 
the frequency enters only in the combination 0^/70, be- 
cause the frequency in g^q appears as e/70 (see Eg. [3T|). 
The equations for the V7-dependent propagators are 



Vi 



oc—eq 



(R; p,Lu) 



■■Veqiuj)Ilioc-eqi^)Veq{uj); (43a) 



Vv7(R; P,t^) = -Veq{uj)Ilvyi^^)Veqiuj). 



(43b) 



We will give a detailed discussion regarding the contri- 
bution of the interaction field to the heat transport in the 
following section. In addition, in Appendix [El we present 
the kinetic equation for an interaction field describing the 
fluctuations in the Cooper channel. 



III. THE ELECTRIC AND HEAT CURRENTS 
IN THE PRESENCE OF A TEMPERATURE 
GRADIENT 

In this section we present the derivation of the heat 
and electric currents in terms of the Green's functions 
as a response to a temperature gradient. We also find 
the heat current induced by an electric field which will 
be used in the next section for the verification of the 
Onsager relations. 

We first consider the electric current. Just like in the 
derivation of the electric conductivity, we start from the 
charge density given by Eq. 1161 and extract the expres- 
sion for the electric current from the continuity equation 
(see Eq. [T7)) . Further calculations using the continuity 
equation requires the equations of motion for the field 
variables which can be obtained from the action pre- 
sented in Eq. [551 At this point, the derivation of the 



electric current as a response to a gravitational field de- 
parts from the one for the electric field, because the dy- 
namics of the field variables is different. According to 
the action given in Eq. [551 the divergence of the electric 
current is 



Vje — 2e lim 

r — 

"V7(r)V - V'7(r')V' 



+ ±y^^r\t'Mr,t)) (44) 



e lim 



2 m 



G(r,t;r',i) 



7(r) / dr,dt,±{r,t;r,,t,)-f{r,)Gir,,t,;r',t) 



+ / dr,dt,G{r,t;r„t,)j{r,)t{r,,t,;r',t)j{r') 



As before, the factor of 2 is due to the summation over 
the spin index. Beside the dependence on 7 through the 
Green's functions and self-energies, the RHS of the above 
equation contains the gravitational field explicitly. We 
may eliminate this explicit dependence on 7 by trans- 
forming to G and E (using the transformation described 
in Eq. 



-G{ri,ti;r',t) 



Vje = -e lim 



2m 

dridti^(r, t; ri , tj_)Q,{r, ,t,;r',t) 
dridtiG{r, t; ri, ti)S(ri, ii; r', t) 



(45) 



In the above equation we used the fact that the product 
7^/^(r)7~^/^(r') in the limit r' ^ r becomes 1. 

We arrived at an equation which is identical to the 
one obtained for the electric current as a response to an 
electric field in Sec. [H Therefore, we get the following 
expression for the current: 



dr'dt' 



v(r,i;r',t')G(r',i';r,i) + h.c. (46) 



The above expression for the current is valid in the 
regime of linear response. One may check that the sec- 
ond non-vanishing term in the expansion described in 
Eq. [201 yields a contribution to the current, je(R) « 
/ dr 1 der 1 (r 1 Vr)2 [^(R; -r 1 , e)^(R; r 1 , e)] < , which is al- 
ready beyond the linear response. 

Next, we shall concentrate on the expression for the 
heat current, j^. Following the prescription used for the 
electric current, we find j/i from the continuity equation 
for the heat density, Q + Vj?i = 0. We start from the 
heat density, which in in the presence of the gravitational 
field is: 



Q(r,t) ==7(r) [h{r,t) - fin{r,t)] , 



(47) 



where the energy and particle densities, /i(r, t) and 
n(r,t), are defined in the absence of the gravitational 
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field. Extracting the Hamiltonian density from the ac- 
tion, we may write the heat density as: 



Q{^,t) = ^ lim 



V7(r)V V'7(r')V' 



2m 



2m 



(48) 



+7(r) (Vrmpir) + 0(r, t)) + 7(r') (T/„„p(r') + 0(r', t')) 

-' S 

-i(7(r)+7(r')) J dr,0(r, i)C/-i (r' ~ r J0(r„ 0^ 



Due to the fact that the average heat density is evaluated 
under the path-integral, Eq. 1481 can be rewritten in the 
following way: 



Q(r, t) — I lim 

r' — *r 









dt'J 



-2S{r-r')d{t-t') 



(49) 



When inserted into the continuity equation, the term in 
the heat density proportional to S{t~t') vanishes. The re- 
sulting continuity equation for the heat current acquires 
a rather simple form: 



Vj/, = -i lim 



fd 


d \ 


{- 








\dt 


dt' / 



(50) 



We have already met the derivative with respect to the 
center of mass time acting on (■0t(r', i')?/;(r, t)) in the 
calculation of the electric current as a response to the 
gravitational field (see Eq. W^i . Following the same route 
as in the transition from Eq.|44]to Eq. l46[ we may express 
the heat current in terms of G: 



) (51) 



X Q{vi,ti;v,t') + h.c. 
de 



2tt 



-dr' 



ev(r,r',e)G(r',r,e) 



+ h.c. 



In the last equality we performed the Fourier transform 
with respect to the relative time coordinate. 

It is worth pointing out that although the interaction 
renormalizcs both the velocity and the Green's function 
in Eq. 1511 for an instantaneous interaction there is no 
direct contribution to the heat current from the inter- 
action propagator. In the scheme developed here, the 
heat density and current were expressed in terms of the 
quasiparticle Green's function alone. This observation is 



consistent with the kinetic equation for V_ that reveals 
that the propagator of the interaction depends on the 
gravitational field only through the quasiparticle Green's 
functions in the self-energy n. The physical picture be- 
hind this two results is connected to the fact that it is the 
quasiparticles that are actually coupled to the heat bath 
and have a well defined temperature, while the interac- 
tion field does not have a temperature of its own. Our 
expression for the heat current is different from the one 
given in Ref. [13 where the heat current has been taken as 
a sum of two contributions, one from the quasiparticles 
and the other from the collective modes. The authors of 
Ref. [13 introduced the two terms because their purpose 
was to construct the kinetic equation in terms of the local 
distribution functions of the quasiparticles and collective 
modes. Since we keep the kinetic equation for the propa- 
gators to be non-local, the quasiparticle Green's function 
entering the current carries the information about the 
collective modes as well. Another difference is related to 
the kinetic equation. In Ref. [l3, the temperature gra- 
dient enters through the derivative with respect to the 
spatial coordinate, rather than the time. Despite all the 
differences, the calculation of the thermal conductivity 
presented in SecfVlproduces the same result as in Ref.hJ. 

Extracting the electric and heat currents from the 
continuity equations is similar to deriving the Ward- 
identities. Therefore, it is not surprising that the ex- 
pressions for the currents contain the renormalized veloc- 
ity. A similar approach for finding the vertex correction 
to the heat current has been applied by Langeril. Our 
expression for the heat current coincides with the first 
term given in Eq. 3.30 of ReffTll. Langer's heat current 
includes also a non-trivial term with a derivative of the 
interaction amplitude with respect to the momentum. In 
the present scheme we succeeded to bypass this compli- 
cation. We believe that in our scheme this term is hidden 
in the V7-dependent self-energy that contains also con- 
tributions in which the gradient of the gravitational field 
enters through the propagator of the interaction. 

For completeness, we derive the expression for the heat 
current as a response to an electric field. In general, the 
procedure is similar to that of finding the heat current 
generated by a gravitational field, but there is one im- 
portant modification. Here, the continuity equation ac- 
quires an additional term due to the work performed by 
the electric field on the electrons: 



Q(r,t) + Vj,,.(r,t) -je(r,t)E(r). 



(52) 



Unlike the heat density in Eq. [171 that is a function of 
the chemical potential /i, in the presence of an electric 
field the heat density is a function of the electro-chemical 
potential. Therefore, the continuity equation for the heat 
density should be written as: 

/i(r, t) - (a* - e^(r))n(r, t) + Vj,,(r, t) ^ -je(r, t) V(^(r) 

(53) 
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(a) 



(b) 



FIG. 3: (a) The heat current vertex for non-interacting elec- 
trons, (b) The heat current vertex in the presence of interac- 
tions is a product of the renormalized velocity multiplied by the 
frequency of the incoming and outgoing legs. 



Using the continuity equation for the n(r, t), one gets: 

h{r, t) - iih{r, t) + V((^(r)je(r, <)) + Vj,,(r, t) = 0, 

(54) 

where the charge current is given in Eq. 1191 The ad- 
ditional term in the heat continuity equation makes the 
expression for the heat current generated by an electric 
field to be gauge invariant. Following the same steps as 
in the derivation of the response to a gravitational field, 
we obtain: 



j,,=i lim + ieip{v) - ^ + «e(^(r') j 



(55) 



r — »r 



dridfiv(r, t; Tj, ti)G(ri, ti; r', t') + h.c. 



Applying the gauge invariant Fourier transform as de- 
fined in Eq.[9l the expression for the heat current in terms 
of the gauge invariant Green's function becomes: 



Jh 



27r 



ei(i-,r',e)G(r',r,e) + h.c. (56) 



Here, the frequency carried by the flow multiplies the 
renormalized velocity. 



We wish to emphasize the most essential feature of the 
expressions for the heat current obtained in our scheme, 
Eqs. [STl and [56l That is, that the frequency factor cor- 
responds to the legs of the renormalized vertex, but not 
to the frequency of the two Green's functions connected 
to the bare velocity inside the vertex, for illustration see 
Fig. El In the diagrammatic technique, one usually starts 
with the bare vertex and then adds the interaction. In 
Appendix [Bl we show that for the simplified Kubo for- 
mula, where the bare heat current vertex is evg, dressing 
this vertex with interactions is not enough to reproduce 
the correct answer. We believe that this is the reason 
why the simplified Kubo formula produces wrong results 
in the presence of interactions. 

We like to point out the similarity in the structure of 
the four currents obtained in our scheme (Eqs. [HI [311 ED 
and I56p . One may notice that after using the proper 
transformations described by Eqs. [51 and [2H1 (for the re- 
sponse to an electric field and a gravitational field, re- 
spectively), all the currents acquire the universal and 
amazingly simple form of Eq. [H In particular, the entire 
dependence of the currents on the external fields is only 
through the renormalized velocity and Green's function. 

From now on we will work only with the tempera- 
ture gradient. This is possible because we already found 
the kinetic equations and the currents and therefore, the 
gravitational field fulfilled its role. We now present the 
expression for the heat current generated by a tempera- 
ture gradient. First we have to adjust the V7-dependent 
parts of the Green's function to describe the response 
to a non-uniform temperature. For that we replace V7 
by VT/T and set 70 = 1 in Eqs. [371 and [21 Next, we 
insert the expressions for Gioc-eq, G'vt as well as the 
VT-dependent velocities into Eq. [51] In the regime of 
linear response, the local equilibrium Green's function 
and the corresponding velocity do not contribute to the 
longitudinal current, because the dependence on the cen- 
ter of mass coordinate makes such contribution vanish 
after averaging over the volume. Thus, the heat current 
is entirely determined by Gvt and Svt: 



J 



Jh 



2T / 27r 



de 



vf{e)^^vf{e)g^(e) - «f (.)3f,(e).f (e) "^^^ 



+ * / ^^^'f (e)4(e) [S^tW(1 - npie)) + ^>A^)nAe)] (4(e) - 4(e)) + , 



a4(e) 



de 



(57) 



We will use this expression to analyze the thermal con- 
ductivity in the following sections. 



IV. ELECTRIC AND THERMAL 
CONDUCTIVITY IN THE FERMI-LIQUID 
THEORY AND THE WIEDEMANN-FRANZ LAW 



In this section we apply the microscopic scheme devel- 
oped in this paper in order to calculate the electric and 
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thermal conductivities in the framework of the Fermi- 
Hquid theory. We demonstrate that the two conductivi- 
ties are related through the Wiedemann-Franz law as it 
should be according to phenomenology. We will consider 
a sufficiently clean system, so that epT ^ 1. On the 
other hand, we will not deal with the hydrodynamic as- 
pect of transport in clean systems and we will assume 
Tt < 1. 

We start with the calculation of the longitudinal elec- 
tric conductivity. For this purpose, we use the expression 
for the electric current presented in Eq. [23l We have to 
perform the averaging over disorder. In the Fermi-liquid 
theory, for a short range disorder, each Green's function 
is averaged separately. Since on average the disordered 
medium is uniform, the equilibrium Green's functions of 
the quasiparticles become diagonal in momentum space 
after averaging: 



93?(Teq(ki?, e 



(58) 



de 



e=0 



From now one, r denotes the renormalized scattering 
time of the quasiparticles. The diagonal form of the 
Green's functions allows us to rewrite the electric cur- 
rent in momentum space. As mentioned above, the sec- 
ond term in Eq. 1231 vanishes for the longitudinal current. 
The last term is proportional to the imaginary part of the 
interaction propagator. This term yields a correction to 
the Drudc conductivity that is smaller— 12£ by a factor of 
(Tr)^. (In fact, this correction is canceled out by other 
terms. However, as we have already mentioned, in this 
paper we will not study contributions that vanish in the 
limit T 0.) Therefore, for the calculation within the 
framework of the Fermi-liquid theory, only the first part 
of Eq. [53] is needed. Let us concentrate on this term: 



J, 



2 J (27r)'^+i 



[(z;f(k,e)+z;f(k,e)) 



de 



(k,6))4(k,e) 



(59) 



-(t;f (k, e) + vtiK e))4(k, e)«f (k, e)4(k, e) 
-<(k, e)(^f (k, e) + «/(k, e))4(k, e)vt{k, e)] 



A non zero contribution in the Fermi-liquid theory is 
generated by products in which one of the two Green's 
functions is retarded, while the other one is advanced 
(see Fig. m. In the above expression, the last two terms 
which contain 5^5^ or 5^5^ vanish after the momen- 
tum integration. Next, we use the fact that the con- 
stant part of the renormalized velocity can be written as 



m IS 



vf(1 -I- 9Kcr^'"^/9^k|^j^^o)' "^^"^re vp = kp/ 
the unrenormalized Fermi velocity. Then, the longitudi 




FIG. 4: The diagrammatic representation of the electric con- 
ductivity in the Fermi-liquid theory. 



nal electric current acquires the form: 

\ 2 „ 



= -2e^E 1 



9a 

dkde 



(60) 



(27r 



\d+l 



dnp(e) „ N 4 /, s 
^^4(k,e)4(k,e) 



where d denotes the dimensionality. At this stage, one 
may neglect the dependence of the Green's functions on 
the frequency in the integration over e. After integration 
over the frequency and momentum, one obtains for the 
electric conductivity: 



'i + a5R(Te,(k,o)/aa|^^^, 

1 - d^aeq {kp,e)/de\ 



m e TIT 



(61) 



IS 

'eg(k 



the 
= kF, 
Thus, 



mass 



m* /m 



renormalized 
we have reproduced the known 



a5RCTeg(k, e = 



where m* 
[1 - aSffe 

expression for the Drude electric conductivity in the 
Fermi-liquid theory, a = e^nr/m*. 

We turn now to the calculation of the thermal conduc- 
tivity. Within the framework of the Fermi-liquid theory 
the calculation of the electric and thermal conductivities 
are practically parallel. Exactly as in the calculation of 
the electric conductivity, the second term in the expres- 
sion for the heat current given in Eq. [57] vanishes. The 
last term is proportional to the imaginary part of the in- 
teraction and it yields a contribution to n/T of the order 
{TtY that can be neglected in the Fermi- liquid theory. 
[We will come back to this term in the next section in the 
discussion following Eq. [5D]] Therefore, the main contri- 
bution to the heat current arises from the first term in 
Eq. 1571 Written in momentum space, this term acquires 
the form: 



= (2^^'^ [(«f (k,.) -f (k,.)) 

x4(k, e){vf{k, e) + «/(k, e))4(k, e) (62) 
-(z;f (k, e) + vt{k, 6))4(k, e)vf{k, e)gf^{k, e) 
-4(k, 6)(«f (k, e) + «/(k, 6))4(k, e)vt{k, e)] . 

Since the last two lines do not contain a pair of retarded 
and advanced Green's functions, they do not contribute 
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to the current. Using the fact that the renormahzed ve- 
locity can be taken outside of the integrals, the longitu- 
dinal heat current can be written as 



VT 

Jh = 2— 



eg 



/ 



dkde 

(2^ 



(63) 



d+l 



de 



(k,e) 



The only difference between Eq. [SO] and the above equa- 
tion is that the heat current contains a factor of while 
the electric current includes the coefficient e^. The ob- 
tained expression is exactly what one expects to get for 
the heat current transported by quasiparticles. It can 
be interpreted as if each quasiparticle contributes to the 
heat current its energy e. This energy is flowing with the 
velocity of the quasiparticle carrying it. 

As a result of integrating over e, the thermal conduc- 
tivity K becomes: 




(64) 



k,0)gfJk,0) 



The remaining integral over the momentum is identical 
to the one encountered while calculating the electric con- 
ductivity. Eventually, we get that the ratio of the ther- 
mal and electric conductivities in the framework of the 
Fermi-liquid, is proportional to the Lorentz number: 



TT 



(65) 



This natural result was obtained in our scheme almost 
automatically, mostly due to the fact that the currents 
have been expressed in terms of the renormalized veloci- 
ties from the very beginning. On the contrary, since the 
Kubo formula starts from the bare vertices, the deriva- 
tion of the thermal conductivity in the Kubo formalism 
is not a trivial taskjii The situation with the simplified 
Kubo formula is even worse. As we show in Appendix iBl 
the simplified Kubo formula generates terms that vio- 
late the Wiedemann-Franz law in the framework of the 
Fermi-liquid. 



V. DIFFUSION CORRECTIONS TO THE 
THERMAL CONDUCTIVITY 

In the previous section we have shown that for a Fermi- 
liquid system the ratio between the thermal and the elec- 
tric conductivities is determined by the Lorentz num- 
ber. One may wonder whether the Wiedemann-Franz 
law still holds for a disordered system when one goes 
beyond the framework of the Fermi-liquid theory in the 
diffusive regime, Tr ^ 1. This question was the sub- 
ject of a long lasting debate. While in Ref. [l^ it was 




FIG. 5: The self-energy before averaging over the disorder. Af- 
ter averaging, it will be decorated with the diffusons. 



concluded that the Wiedemann-Franz law remains valid, 
later studies show that, in fact, it is violatedJ^ii^iiii^^ 
To find the answer to this question using the scheme pre- 
sented here, we study the leading order corrections with 
respect to {eft)'^ to the thermal conductivity. In the 
electric conductivity these are the Altshuler-Aronov cor- 
rections!^ In two dimensions, the Altshuler-Aronov cor- 
rections logarithmically diverge as the temperature goes 
to zero, Sa oc e^lnl/Tr. Correspondingly, we wish to 
examine the singular corrections to the thermal conduc- 
tivity for which k/T logarithmically diverges in the limit 
T ^ 0. We will show that some of the corrections, origi- 
nating from the long range Coulomb interaction, violate 
the Wiedemann-Franz law. These corrections emerge 
from the third term in Eq. [57l which is related to the 
imaginary part of the polarization operator. While in the 
Fermi- liquid calculation they are proportional to (Tr)^, 
when dressed by the diffusion propagators (diffusons), 
they are dramatically increased and do not contain this 
smallness anymore. Our result agrees with the one given 

in Refs. 

To study the corrections to the transport coefficients in 
the diffusive limit, Tt <^ 1, one has to perform properly 
the averaging over the disorder. For that one has to 
construct the diffusion propagators in the following way: 

(<?f/(r, n, e)4^«(r„ r, e - u;)),mp (66) 

= I rfr'i?«'-4(r,r',..)5f/(r',ri,e)4«(r„r',6-c.), 



where the diffusion propagator is: 



{67) 



Here, D is the diffusion coefficient. As we shall see, we 
need to consider contributions which may include up to 
four diffusons. To obtain the lowest order corrections in 
(eft)^^, all the diffusons must carry the same momen- 
tum. This is because each integration over the momen- 
tum of the diffusion propagators produces a small factor 
{eft)^^. This argument implies that all the diffusons 
must be affiliated with the same self-energy. Therefore, 
we generate the singular corrections by choosing one of 
the self-energies to be decorated by the diffusons. When 
the chosen self-energy is the equilibrium one, we denote 
it by a'^'^ff , while for the VT-dependent self-energy we 
use the notation fl'^'lj/ . 
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Next, we expand the expression for the heat current 
given in Eq. [ST] with respect to the choosen self-energies. 
In the first term of the heat current, we expand both the 
equihbrium Green's functions and velocities with respect 
to a'^^ff. For example, the expansion of the Green's func- 
tion yields geq = 9 + ga-'^^^-^g, where g is the equilibrium 
Green's function which already incorporates the Fermi- 
liquid renormalizations. As has been mentioned before, 
the second term of the heat current does not contribute 
to the longitudinal conductivity. Unlike the calculation 
in the previous section, the third term in Eq. [57] is no 
longer negligible. In order to compensate the smallness 
initially associated with the third term, it is Evt that 
must be decorated by the diffusons. Since the number 
of integration over the momentum of the diffusons is re- 
stricted to one, the self-energy chosen to contain the dif- 
fusons can have only one effective interaction amplitude 
as illustrated in Fig. [5] [Note, that the expansion is with 
respect to (ej-r)"^, and not over the interaction.] For the 
Coulomb interaction, the effective interaction amplitude 
can be approximated by the propagator V. Hence, 

E<'>(R;p,r) =zG<^>(R;p,T)F<'>(R;p,r); (68) 
E^///R; p, r) = jG«^^(R; p, r)V<{IL; p, r) 

+ iG>(R; p, r)l/^^^(R; p, r). 

The interaction propagator is described by the infinite 
geometrical series presented in Fig. [6] Each term in this 
series is a convolution in space. Symbolically, the result 
can be written as: 



V{uj)= U-^ +IV{uj) 



(69) 



The polarization operator 11 has the following analytic 
structure: 



n 



■ (R; p, r) = zG<'> (R; p, r)G>-< (R; -p, -r); (70) 



R,A/ 



(R; p, r) = zG«'^(R; p, t)G< (R; -p, -r) 



+ iG<(R;p,r)G^'«(R;-p,-T). 

In Eq. [68l S denotes either the equilibrium or the VT- 
dependent self-energy. For Evt we should consider all 
possibilities to linearize the expression in Eq. [68] with 
respect to VT including the Green's functions inside the 
polarization operator. 

We are now fully equipped to study the leading order 
corrections to the thermal conductivity induced by dis- 
order. We separate the corrections to the current into 





(a) 



(b) 




(c) 



FIG. 7: The terms with a horizontal interaction line after aver- 
aging over the disorder. 



three groups. The first group includes all the contribu- 
tions which arise from the expansion of the equilibrium 
Green's function in the first term in Eq. 1571 



h. 



de 

dT J — 



w^(r6,ri) 



X 9 (ri,r2,e)CTdi//(r2,r3,e)5 '^{v^,Y^,e)v' {v^,v.^) 



b (i-5,r6,e) -5^(r5,r6,e)] +c 



(71) 



Here, we used the fact that for the longitudinal currents 
[i = j) the contribution from the term v"^ g^a^^^ ^g^v"^ gn 
is the same as the one from g^v^ g^a^^^ ^g^ . After per- 
forming the average over the disorder in j/j /jo^: we get: 



3h ho 



2i- 



.V,T 



T 



{npie) - npie - w)) 

2 
f 

d 



dnp{—ui) dnp{€] 



np(— a;) 



duj de 
g:(5^(k, e))2g^(k -q,e-Lu) [g«(k, e) - g^ik, e)] 



2TlVTd 



dk 



^(g^(k,e))V(k,e)(g^(k',e))= 



xg'^(k' - q, e - w) + 



2TivTd J {2ttY m 



'^■'^(5^(k,e))^ 



x/(k - q, e - ^)(g«(k', e)fg^{k! - q, 6 - u:)Dr{ci, u\ 



+ c.c. 



(72) 




FIG. 6: The series describing the interaction propagator. 



The contributions to the current hor can be inter- 
preted in terms of the diagrams with a horizontal inter- 
action as shown in Fig. [7] The only difference between 
the expression written above and the corresponding con- 
tributions to the electric conductivity is that should 
be replaced by — e^. This is the proper place to remind 
that our scheme does not require a diagrammatic calcu- 
lation, rather all the contributions are generated using 



16 




(a) 



(b) 



FIG. 8: The diagrams with a vertical interaction line. The 
obvious counterpart of diagram (b) is not shown. 



the quantum kinetic equation. In this method the ana- 
lytic structure of each of the terms and their numerical 
coefRcients are determined by the kinetic equation. We 
give a diagrammatic interpretation of the different term 
for the purpose of illustration only. 



In the second group we collect terms related to the 
renormalization of the velocity. These contributions orig- 
inate from the first and the third terms in Eq. [ST] In the 
first term we expand one of the velocities with respect to 
In the last term, we consider the contributions in 
which the the temperature gradient enters S^^''^ through 
the quasiparticle Green's function. [We relate the other 
contributions to the Coulomb drag that will be discussed 
afterward.] Performing these operations, we get: 



Jh ve 



dediv 
(2^ 



dr 1 . . . drieu* (r4 , r 1 (r 1 , ra , e) 



X [(1 - nF{e))V<{r^,r,„uj)G^j.{r2,r^,e - uj) 
+nF (e) V> (r^ , , u)G^j, (r^ , T;, , e - tj)] 

X [.g^(r3,r4,e) -5'^(r3,r4,e)] j ^dr,...dr^e^ 

X [ASvfir^, ri, e)g^(ri, r^, e)w^(r2, r3)5'^(r3, e) 
-3(5wf (r4, ri, e).9^(ri, r^, e)w^(r2, r3)g^(r3, r^, e) 

+ c.c. (73) 



Here, 5v^'^{r,r' ,e) = -i{n ~ r'^)a^^ff{r,r' , e). Since 
in this calculation Gvt includes only the Fermi-liquid 
renormalizations, it is described by the last term 
in Eq. Ei GvT(e) = -ieVT/T[dg{e)/dew{e)g{e) - 
g{e)v{e)dg{e) / de]. As a result of averaging over the dis- 



order we get: 

ViT f dkde dadw ,^d, 

dnp{~uj) 



Jh 



I ver 

4eV,^(q,c.) 
dnp{e 



de 



np(—LL}) 



+ ecu [ni?(e) — nF{e — uj)] 



dnp{uj) 



X [y«(q,^)-F4(q,c.)] 



x(g-4(k-q,e-u;))2 + 



2tivt J {2'kY w? 
X5^(k - q, e - ^)(5^(k', e - uo)fg''{\^ + q, e)Dp{ci, oj] 
+ c.c. 

In the transition between the last two equations as well 
as between Ea . [73l and Eq. [74]we used the standard iden- 
tities for the product of Fermi distribution functions: 



dnF{e) 
de 



-npU 



dnp{—u)) 
duj 

np{—LL}) 



[nF{e) -npie- uj)] 



dnF{e) 
de 



(75a) 



nF(e){\ — npie — u)) = np{uj)[nFie — u>) — nF{e)]. 

(75b) 

The diagrams corresponding to the second group are 
presented in Fig. [5] One can observed that the second 
group contain the contributions with a vertical interac- 
tion line. The terms containing in Eq. [74] have their 
counterparts in the corrections to the electric conductiv- 
ity. Together with the first group they give a correction to 
the thermal conductivity which satisfies the Wiedemann- 
Franz law: 



dqdio ^dnpie) d{ijjnp{~uj)) 
-dee 



(27r)'*+i 



de 



Dq'Vj}{q,u;)Dl,{q,u;y 



c.c. = —SaAA, 



(76) 



where Sctaa is the Altshuler-Aronov corrections to the 
electric conductivity. The additional terms that are pro- 
portional to ecu in Eq. [74] are responsible for the devi- 
ation from the Wiedemann-Franz law. Only the dia- 
gram with two diffusons gives a singular contribution 
which in two dimensions is accumulated in the region 
of small momenta: uj/ Dugcreen < Q < \/^/D, where 
1/Kscreen = l/(27re^z^) is the inverse screening length in 
d ~ 2. Eventually, the non Wiedemann-Franz law cor- 
rection to the thermal conductivity is 

X ■ r, ^ I ^^^^ 

dnp{bj) 

doj 
12 



(27r)3''^"'^ ['^F(e) - npie- uj)] 



K^(q, uj) - F4(q, io)] [i?|j(q, uj) + D\{q, cu)] 



T 



(77) 
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Although both Skwf and Sunon-WF logarithmically di- 
verge in d = 2, the origin of the singularities are differ- 
ent. The logarithmic correction that do not violate the 
Wiedemann-Franz law accumulate over a wide region of 
momenta y^T/D < q < 1/vpT. 

The third group corresponds to the Coulomb drag, see 
Fig. [HI These terms are generated when the temperature 
gradient enters S^^''^ through the interaction propaga- 
tor. To exploit the symmetry related to the Coulomb 
drag, we shall use the fact that the interaction field cj) 
is real. Therefore, the corresponding propagator satis- 
fies the relation: V{r,t;r\t') = -i{Tc{(f>{r,t)(j){r' ,t')}) ^ 
V'^{r',t';r,t), and we may write the lesser and greater 
components of the interaction propagator as: 

V<-> (p, r) = 1 [V<^> (p, r) + F>'< (-P, -r)] . (78) 

Consequently, the self-energy given in Eq.[5S]can be split 
into two parts: 



2 J 2tt 

+G<'> (R; p, e + uj)V>^< (R; -p, uj) 



(79) 



This way of writing the propagator V and the self-energy 
turns out to be highly useful in the derivation of the 
Coulomb drag terms as well as for the proof of the On- 
sager relation presented in the next section. 

A full analysis of the Coulomb drag contributions to 
the thermal conductivity, starting from Eq. 1791 is pre- 
sented in Appendix. [C] We show there that in the ther- 
mal conductivity, unlike the electric conductivity^, the 
Coulomb drag can be decorated by four diffusons (see 
Eg. IC9p . After averaging over the impurities, the correc- 
tion to the heat current from the Coulomb drag is 



Jh drag 



ViT f dqduj 



d f dkde 



^—^\Ve,i<l,u;)\ a;- 



, dnp{ui) 
doj 



(80) 



dq,J (27r)' 



+T["-F(e - w) - riFie)] [g^(k -q,e-uj) 



X5^(k, e)D^iq, to) - 5^(k - q, e - c.)g«(k, e)D'\q, u;)\ 

Owing to the structure of this term, there are no diver- 
gencies related to the region of small momenta indicated 
above. Moreover, the integration over the frequency uj 
is restricted to \uj\ < T. Therefore, this contribution to 
k/T is regular. [Note, that the structure of the Coulomb 
drag term presented here differs from the one in Ref. [28l 
obtained from the simplified Kubo formula.] In the pre- 
vious section we argued that in the framework of the 
Fermi-liquid theory the third term in Eq. [57] generates 
contributions that are proportional to (Tr)^. Using the 
expressions in Eqs. [771 and [501 one may check that in the 
absence of the diffusons the contributions related to the 
third term indeed acquire this small factor. 
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FIG. 9: The drag diagrams after the averaging over the disorder. 



To conclude the section, we return to the singular 
correction to the Wiedemann-Franz law presented in 
Eq. 1771 The infrared divergency of this kind does not 
occur in the electric conductivity due to gauge invari- 



ance 



To understand this argument, one has to 
go back to Eq. [Sj describing the action of interacting elec- 
trons. Writing the action in momentum space and con- 
sidering the limit q — > 0, one may neglect the depen- 
dence of the quasiparticle Green's functions on the mo- 
mentum transferred by the interaction. Then, after the 
integration over the transferred momentum, the field 4> 
that describes the effective electron-electron interaction 
becomes only a function of time. An interaction field 
of the form (j){t) can be gauged out by redefining the 
quasiparticle field ip{r,t) e'' J'' '^*''''^-*">-ip{r,t). In the 
presence of a temperature gradient the situation is more 
complicated. The correction to thermal conductivity vio- 
lating the Wiedemann-Franz law are proportional to the 
imaginary part of the interaction propagator V. In other 
words, these terms are sensitive to the decay of the inter- 
action into particle-hole pairs. Therefore, once the inner 
structure of the interaction becomes important, gauging 
out the interaction field on the level of the action is no 
longer justified. 



VI. ONSAGER RELATIONS 



In Sec. IIIIl we derived the expressions for the electric 
current as a response to VT and the heat current gener- 
ated by E. The two expressions were obtained indepen- 
dently from each other. On the other hand, these two 
currents must be related through the Onsager relations. 
In this section, as an additional test for our scheme, we 
verify that the expressions given in Eas. l46l and l56l indeed 
satisfy the Onsager relation. 

The Onsager relations22, connect between the off- 
diagonal elements of the conductivity tensor: cTy (B) — 
(7J^{-B), Ky(B) = Kj^(-B), and 5jj(B) = Taj.(-B). 
In the absence of magnetic field, the Onager relations 
reduce to ctxx = Ta^x- Here, we restrict our demonstra- 
tion of this relation to a system of electrons interacting 
through the Coulomb interaction as has been discussed in 
the previous sections. [In this paper, our treatment of the 
electric and heat currents does not include the modifica- 
tions needed to account for the effect of a magnetic field. 
The generalization of the present scheme for calculating 
the thermoelectric transport coefficients in the presence 
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of a magnetic field is straightforward, and was applied by 
us in the analysis of the Nernst effect in disordered films 
above the superconducting transition temperature^^.] 

To find the two thermoelectric currents, we have (i) 
to insert the E-dependent velocity and Green's function 
into the expression for the heat current given in Eq. 1561 
and similarly (ii) to insert the VT-dependent velocity 
and Green's function into the electric current described 
by Eq. |46l [As has been already argued in the deriva- 



tion of the thermal conductivity (see the discussion above 
Eq. I57p. the local equilibrium Green's function does not 
contribute to the longitudinal currents.] Due to the sim- 
ilarity between Ge and Gvt and the common structure 
of the four currents (Eqs. [21 SI EH and [SS]) , aU the re- 
sponses determining the longitudinal components of the 
conductivity tensor can be written as: 



-<(e)^;f (e)4(6)^;/(6)] + ^ J ^^^^^{e)vfie)gf^{e) [E<(6)(1 - np{e)) + Y.>{e)np{e)] (e) - 4(6)) -h c.c. (81) 



Here, F(e) is equal to eE for the response to an elec- 
tric field and to eVT/T for the response to a temper- 
ature gradient. For the electric current, Xe(e) = 
while for the heat current, Xh{^) ~ £• Comparing the ex- 
pressions for je(Vr) and jh(E), one may immediately 
see that the contributions described by the first term 
in the above expressions satisfy the Onsager relation, 
jl/Ei — —Tjl/^iT. Therefore, it remains to show that 
the same holds for the last term. The dependence of the 
last term in Eq.[8T]on the external perturbation F enters 
through the self-energy. As has been mentioned in the 
previous section, we need to consider all the possibilities 
to take one of the propagators (G or V) in Ep to depend 
on F. 

In general, the fulfilment of the Onsager relation de- 
mands microscopic reversibility, which in our case implies 
that G(r,r',e) = G(r',r,e) and V{r,r',e) = V{r',r,e). 
The Onsager relation is satisfied if by reading the contri- 
butions to je(Vr) described in Eq. W\\ from right to left 
instead of left to right, one gets jh(E), and vice versa. 
Besides the microscopic reversibility, in order to get the 
desired relation the thermoelectric currents should have a 
proper symmetry related to the frequency. There should 
be a well defined symmetry (for each of the currents) un- 
der the exchange of the frequencies carried by the Green's 
functions attached to the two velocity vertices. The sec- 
ond velocity vertex appears in the last term of Eq. 1811 
because, as can be seen from Eqs.[T51and 15^1 the external 
perturbations are accompanied by the velocity. 

We will show now that this symmetry is embedded in 
the structure of the combination Z entering the last term 
in Eq. [ST] 



Z = S<(e)(l - npie)) + Y.>{e)np{e) 



(82) 



Recall that Ep contains one F-dependent propagator, Gp 
or Vp, that by itself can include a self-energy that de- 
pends on the external perturbation, Sp or Hp, i.e., Ep 
is determined iteratively. We start our analysis of Z at 



the point when the iterative process was already termi- 
nated. In other word, we will use the expression for Sp 
at a stage when the Green's function depending on the 
external perturbation is equal to: 



(83) 



de 



v(e + f7)geg(e + f^) 



-5eg(e + f^)v(e + f^) 



de 



Here the argument of the Green's function reflects the 
fact that in the proof of the Onsager relation we will 
apply the standard trick used in the derivation of the 
Ward identities. Namely, we will arrange the arguments 
of the Green's functions inside Sp (e) in such a way that 
they all include e (even those inside the loops). 

Any given contribution to S^'^ is a sum of few terms 
with a different analytic structure of the propagators. 
Each of these terms can be cut into two pieces in such 
a way that the propagators along the cut are the lesser 
or greater components only. [For a self-energy with two 
crossed interaction lines all possible cuts are shown in 
Fig. [TOl] Therefore, each term can be written as a prod- 




la) 



(b) 



(c) 



FIG. 10: The self-energy in the second order expansion with 
respect to the interaction. 
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uct A- B, where B includes all the propagators along the 
cut, while all the rest of the propagators are collected 
in A. Then, each term in — A ■ has a counter- 
part in = A ■ B^ with the same A but "opposite" B. 
"Opposite" means that if we substitute in B^ the lesser 
propagators by the greater ones and simultaneously all 
the greater by lesser, we get B^ . 

To demonstrate how this prescription works, let us look 
at the expression corresponding to the diagram presented 
in Fig.[ini;b): 

S<'>(r,r',e) = -i J ^cir,dr.G«(r,r„.-.) 

X G<->(ri, r„ e~uj~ Lu')V<^>{r, r„ c^)V^<'> (r„ r', lo') 

X G'^(r„r',e-^')- ^ J ^|^dridr,G^(r, r^, e + w) 

X G<^> (ri, r^, e + uj + Lu')V>'<{r2,r, lu)V>'< (r', r,,Lu') 
xG'^{r2,r\e + uj'). (84) 

In this example A = G^G^, B< = G<V<V<, and 
= G-^V^V^ . The above expression contains two 
contributions in which the frequencies uj and lo' appear 
with opposite signs. To get the two contributions we 
split the lesser and greater components of the interac- 
tion propagator into two pieces according to Eq. (75] We 
will use this representation of V throughout the proof of 
the Onsager relation in order to write the self-energy as 
a sum of two terms in which the F-dependent Green's 
function appears either as Gf{£ ~ ^) or GF(e + fl). 

Let us return to the discussion of the general properties 
of Z. After a pair, AB^ and AB^ , is inserted into Eq.[5^ 
one has to linearize it with respect to the external pertur- 
bation. Recall that in the end of the iterative process, F 
enters the self-energy through the quasiparticle Green's 
functions Gfs- When this Green's function belongs to 
A, the corresponding contribution to Eq. [82l ^f[(1 — 
np{e))B^g{e) + np{e)B^g{e)], vanishes. This is because 
the total frequency transferred along the cut is equal 
to e, and at equilibrium B(e) has the same properties 
as any other fermionic propagator in the Keldysh tech- 
nique. Therefore, the only non-vanishing contribution is 
of the form Z = ^eg[(l - nF{e))B< + nF{e)B>]. Be- 
fore we start analyzing Bp, we wish to note that we may 
extend A^q to include all the terms with the same cut. 
The Green's function through which the external pertur- 
bation enters Bp, GFs(e ± ^1), is described by Eq. [531 
The rest of the propagators in B which are at equilib- 
rium carry altogether the bosonic frequency and will 
be denote as Peq{^)- Then, we can write B as -B^(e) = 
G<, (e - n)Peq {n)np (Q) + G<^ (e + n)Peg {n){i + np {n)), 
and correspondingly, B^{e) = Gp^{e — il.)Peq{^){l + 
np(r2)) -I- Gpj,(e + Q)Peq{^)np{Q). In the last two iden- 
tities we assumed that the F-dependent Green's function 
in B^ is Gp^, while in B^ it is Gp^. [The possibility 
that the external perturbation enters B^ {B^) through 
Gpj, (Gpg) will be discussed below.] Inserting the above 



expression for into Z, we get: 

Z = [{1 - np{e))np{n)G<^{e - n) (85) 
+ np{e){l + np{n))G>^{e - n)]A,qP,q{n) 
+ [(I - np(e))(l + np{n))G<^ie + n) 
+ np{e)72p{n)G>^{€ + n)]AeqPeqin). 

Applying the Keldysh rules on the product of matrices 
in Eq. [83l we obtain that G^^(e - fl) ^ np{t - n)X + 
dnp{e-n)/dey, and G^^{e-n) ^ {np{e-n) - 1)X + 
dnp{e — ^)/dey. From the identity given in Eq. [tSH it 
follows that the terms proportional to X in Eq. [85] vanish. 
Then, using the identity for the product of distribution 
functions given in Eq. I75al we get that the discussed 
contributions to the thermoelectric currents arc of the 
form: 

,88, 

X [/e,ft(e)/iF(e - r!) + /e.,,(e - r!)/2F(e)] . 

In the second term we shifted the frequency e by fi. This 
shift of the frequency affects I2 in such a way that under 
the condition of microscopic reversibility, the functions 
Ii and I2 transform one into another when read in the 
opposite directions. Here we rely on the fact that all the 
contributions with the same cut have been included in 
Ii and 12- Then, it follows from Eq. [SSjthat all the con- 
tributions to the thermoelectric current under discussion 
satisfy the Onsager relation. 

It remains to examine the case when the external per- 
turbation enters B"^ through Gp^ and B^ through Gp^. 
Then, B can be written as B< = G^^(e - n)Peq{n){l + 
np{n)) + Gp^{e + n)Peq{^)np{il.), and correspondingly, 
B> = G<,(e - n)Peq{n)np{n) + G>^(e + n)p,q{n){i + 
np{Q)). One can check, in the same way as before, that 
the contributions to the current from such terms are also 
described by Eq. [86l Thus, we have shown that the last 
term in Eq. [HH is compatible with the Onsager relation. 
In Appendix [D] we give representative examples illus- 
trating the general arguments presented here. 

To summarize, we proved that the expressions for the 
longitudinal thermoelectric currents obtained using the 
quantum kinetic approach (presented in Eq. I8ip satisfy 
the Onsager relations. We use the Onsager relations as a 
check for the correctness of our scheme. In the presence 
of a magnetic field the transverse thermoelectric currents 
contain additional contributions and, hence, the proof of 
the Onsager relations should be modified accordingly. A 
demonstration of the validity of the Onsager relations for 
the thermoelectric transport coefficients in the presence 
of a magnetic field for electrons interacting in the Cooper 
channel will be given in a separated manuscript. 

VII. SUMMARY 

We developed a new scheme for studying thermal and 
thermoelectric transport in interacting electron systems 
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in the presence of disorder. The kinetic equation in the 
presence of a temperature gradient was derived directly 
from the action. One of the novel aspects of the scheme 
regards the expressions for the four currents correspond- 
ing to the different components of the conductivity ten- 
sor, see Eqs. \22 \ l46 l ISD and [56l These currents, derived 
from the continuity equations for the charge and the en- 
ergy, share a uniform and compact structure summarized 
in Eq. [4l This equation reveals that the expressions for 
both the electric and heat currents include the renormal- 
ized velocity. The frequency factor in the expression for 
heat current corresponds to the legs of the renormalized 
velocity as illustrated in Fig. [3] (but not to the frequency 
of the two Green's functions connected to the bare ve- 
locity inside the vertex). This observation, which is the 
main advantage of the scheme presented here, implies 
that in the heat current the flow of energy occurs with the 
renormalized velocity. As we demonstrated in Sec. IIVI 
this structure of the heat current guarantees that the 
Wiedemann-Franz law is satisfied for Fcrmi-liquid sys- 
tems. 

In this paper we considered the Coulomb interaction, 
which is instantaneous in time. Applying our scheme for 
electron-electron interaction mediated by phonons that 
by themselves carry energy is straitforward. An example 
for electron-electron interaction of a different kind has 
been briefly described in Appendix [E] There we con- 
sider an interaction mediated by superconducting fluc- 
tuations. In this case, the electric current acquires an 
additional contribution because the fluctuations in the 
Cooper channel carry charge. Nevertheless, this addi- 



tional contribution does not ruin the general structure of 
the current. 

The main strength of our scheme is in its generality 
and compactness. This is the reason why we believe that 
it is an adequate alternative to the Kubo formula, which 
for the thermal transport is rather cumbersome. 
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APPENDIX A: THE ELECTRIC FIELD 
DEPENDENT GREEN'S FUNCTION 

In this appendix we wish to present additional details 
of the derivation of the quantum kinetic equation in the 
presence of an electric field. In particular, we show how 
to obtain the electric field dependent Green's function 
which solves the kinetic equation, sec Eq. [121 

We start from the Dyson equation for the Green's func- 
tion of the quasiparticles given in Eq. [T] After perform- 
ing the gauge invariant Fourier transform introduced in 
Eq. [HI we obtain the quantum kinetic equation for the 
gauge invariant Green's function: 



i d 
2df 



2m dp^ 



eE 92 



8™ d'R? 8m de^ 



Am dedR 2m dpdH 2m dpde 
X G(R, T; p, e) = 5{p) + j dt,dr,dTde,de,e"'''-"'''-^-''^-"''''± (R + r,/2, T + tj2; p-r„e. 



V,rr,p{-R + p/2) 



epE' 



(Al) 



The proof that the Green's function transformed in this way is gauge invariant as well as other useful technical details 
can be found in Ch. 7 of Ref. [U The quantum kinetic equation is considerably simplified if to restrict the calculation 
to the steady state response (i.e., a DC electric field). Then, the dependence of the Green's functions and self-energies 
on time is only through the relative time-coordinate r. After expanding the phases in the RHS with respect to the 
electric field, the kinetic equation acquires the form presented in Eq. 1101 

Next, we write the Green's function as G = geq + Ge and, similarly, we replace the self-energy by S = deq + Se- 
As the properties of the equilibrium Green's function have been already described in Eq. [T2j here we shall examine 
only the electric field dependent part of the Green's function. For this purpose, we collect all the terms in Eq. \10\ 
which are linear in the electric field: 



2m 



e + h /i - Vimp Ge(p, e; imp) - / drideqip - ri, e; imp)GE(ri, e; imp) 



(A2) 



'5(p-rO 



d&eq (p - ri,e;imp) 



(5(p-ri) d 



X 9eq (fi, e; imp) + / drjl^ {p 



de 

ri,e;imp) g^q {r„e;imp) 



m 



Or I 



+ (p - ri)(Teq (p - ri, e; imp) 



d_ 

de 
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As we have already mentioned, the entire dependence of the gauge invariant Green's functions and self-energies 
on the center of mass coordinates is due to the scattering by the impurity potential. This dependence has been 
incorporated in our notation into imp. We will use the fact that the equilibrium Green's function is g^q{p, e; imp) = 

e + ^ + M ^ Virnp — <^eq{Pi ^] imp) (scc Eq. [TT|) . Then, we may rewrite Eg. IA2l as: 



d^i9eq{p - e; ^mp)G^{v^,e\ imp) = / dv^dt^ [p - Ti, e; imp) g^g {r„e; imp) 



S{p-r,) - 



d&eq {p - fi, e; imp) 



de 



Sjp-r,) d 
m dvi 



(p - ri)(7eg (p - ri, e; imp) 



(A3) 



One may check that the following identities for the equilibrium Green's function hold: 



dgeqip,£;imp) 
de 



dv^dv^geq {Yi,c,imp) 



5{p - ri - ra) 



da^q (p-r^ -~r2,e;imp) 



de 



geqir2,e;imp); (A4) 



-ipgeq{p,e;imp) = ~i J dridr2geq (ri,e;imp) 
= 5(e)v(e).g(e)- 



S{p - ri - r^) d 
m dvn 



(p - ri - r2)^eq (p - ri - r^, e; imp) 



geq {r2,e;imp) 



Using the above identities, one may convert Eq. IA3l into 
Eq. [13] describing Ge- 

The expression for Ge is used in order to transform 
from Eq. [51] to the final formula for the electric current 
presented in Eq. [23l Before inserting Eq. [13] into Eq. [22l 
one should recall that the rcnormalizcd velocity which 
appears in the expression for the current also depends 
on the electric field through the self-energy term. There- 
fore, the current linearized with respect to E contains 
two terms: 



ze 
' 2" 



dr'dt' 



Ve,(r,t;r',i')GE(r',i';r, t) 



+VE(r, t; r', t')geq{r', t'; r, t)]^ + h.c, 



(A5) 



where VE(r, t; r', t') = — i(r — r')SE(r, t; r', t'). One may 
notice that in the limit r' ^ r and t' ^ t (i.e., under the 
trace) the term with the electric field dependent velocity 
can be arranged as follows: 

- i{r ~ r')SE(r, i; r', t')g,q(y' , t'; r, t) (A6) 
= SE(r, t; r', t')i{r' - r)geq{r' , t'; r, t) = -^Egeq-Vcqgeq- 



Finally, one should replace Ge in Ea. lA5l with its explicit 
expression given in Eq. [13] To get Eq. [23l it remains to 
extract the lesser component. For this purpose, we use 
the fact that the current is diagonal in the energy (con- 
volution in time), and thus in the basis of the retarded, 
advanced and Keldysh components the product of matri- 
ces obeys the following rules: 



APPENDIX B: CALCULATION OF THE 
THERMAL CONDUCTIVITY USING THE 
SIMPLIFIED KUBO FORMULA 

In this appendix we wish to examine the thermal con- 
ductivity calculated using the simplified Kubo formula 
presented in Eq. [2] Let us consider one representative 
term in the perturbation expansion of the Kubo formula 
with respect to the electron-electron interaction. For 
this purpose, we choose the second order term drawn 
in Fig. 111! In the Fermi-liquid theory, the main contri- 
bution to the transport coefficients arises from diagrams 
which contain one pair of Green's's functions with the 
same arguments but opposite analytic structure, i.e., one 
is retarded while the other is advanced. Then, the quasi- 
particles corresponding to these two Green's's functions 
are both on the mass shell. This pair yields a contribu- 
tion to the transport coefficients that is proportional to 

T. 

According to Eq. [2l the expression corresponding to 
the diagram presented in Fig. II II is 



K = -—T^ y y 



ild ^-^ ^-^ m \ 2 



X 1 I e„ + 



^ ) 5(k - q, e„ - lj„j + Vl) 



X g(k - q, e„ - lo„x)V(s\-, c^m)5(k, e„ -|- ri)g(k, e„) 

X !/(q',w^)5(k-q',e„ - w^/ f^)5(k - q', e„ -w™'). 

(Bl) 
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(A7) 
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= A^B^; 




{A 


B 




^ A"" B'-" + A'' B^ . 





In addition, we use the relation A< = {A'^ ~ A^- + A^) /2. 



Notice, that in the above expression we use the heat cur- 
rent operators of free electrons. Here, e„ — 2nT{n + 1/2) 
and ojm = 2iTTm are the Matsubara frequencies, and 
the Green's function of the quasiparticles is ^(k, e„) = 
[ie„ — + ^sign{en)]~^ ■ In search for the DC conduc- 
tivity, we should expand the expression to a linear order 
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in the frequency of the external field, fl. Keeping only the terms with one g^g^ pair (R-A section), we get: 
I 



k,q,q' -n<e<0 

h - q, 



oo oc 



E E+E E+E E+EE 



— — oo uj„,/ — — oo cj„ 



-OcJ„,/— — oo UJm—^'iL 



(B2) 



e„ - UJ,n + ^ j - Wm' + ^ ) V^(q, CJ„i)y (q' , LJ„i/ ) 



/fci / ri\ fc, -9' 
+2— e„ + — I e„ - w,„' 

TO \ 2 / TO 



— ) T/^(-q, -uj„i)V{c[ + q,ujm+ t^m') 



X 5(k - q, e„ - + ri)g(k - q, e„ - w„)g(k - q', e 

r 



g(k, e„ + l^)5(k, e„) 
ri)g(k-q',e„ - Wm')- 



In the above expression, the first line describes all the 
cases where the R-A section is in the middle of the dia- 
gram, while in the second line the R-A section is located 
in one of the sides. The arguments in Eq. IB2I are ar- 
ranged in such a way that each time the frequency of the 
effective R-A section is denoted by e„. Then, the sum 
over the frequency e„ is restricted to a narrow window of 
the width fl. Therefore, in the limit ^ 0, the external 
frequency in the Green's functions, the heat current ver- 
tices and the sum over the frequencies and tOm' can 
be set to zero. 

In the next step we perform the analytic continua- 
tion from the Matsubara frequencies to the real axis, 
replacing the sums over the frequencies with integrals. 
In the present calculation, the only difference between 
the thermal and the electric conductivities is that in- 
stead of the frequencies in the heat current vertices the 
electric current vertices give a factor of — e^. There- 
fore, for the Wiedemann-Franz law to be valid, the 
integral over the frequency e should be of the form 
/ c?e(9tanh(e/2r)/9e)e^. We are going to check whether 
all the contributions to the thermal conductivity given in 
Eq. IB2I are indeed of this form. With this in mind, we 
separate the contributions to the thermal conductivity 
into two group. The first contains all the terms in which 
the R-A section is connected to a bare vertex. Let us 
concentrate on the e integration, which is the same for 
all the terms in this group. In the limit Tr 1, we may 
neglect the dependence of all the Green's functions on e. 



Then, the integral over e reduces to: 



, 9tanh(e/2T) , 
de — e (e - ) . 



(B3) 



Only the part in the above integral results in a finite 
contribution, while the rest being an odd function of the 
frequency vanishes. Naturally, the contribution of this 
group acquires the following form: 



E 



de atanh(e/2T) h 



de 



e^gn{k,e)gA{h,e) (B4) 



d^a^ Ck, e) 

^ — It — ^ 
oh 

where 172 is one of the contributions to the self-energy 
with two interaction lines (see Fig. [TWb)). Clearly, such 
a contribution is consistent with the Fermi-liquid expres- 
sion for the thermal conductivity presented in Eq. 1631 
and as such satisfies the Wiedemann-Franz law. 

We turn to the second group that, as we will see, ap- 
pear to be inconsistent with the Wiedemann-Franz law. 
Now the R-A section is located between the two inter- 
action lines. Expanding with respect to the external 
frequency and performing the analytic continuation, one 
gets: 



li,E,+£l 






n)„.+ Sill) 

(a) (b) 

FIG. 12: The contributions to the self-energy of the quasipar- 
FIG. 11: A contribution to the thermal conductivity calculated tides relevant for the discussion of the diagram presented in 
using the simplified Kubo formula given in Eq. [S] Fig. 1111 
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1 



k,q,q' 



duj' , 

cotn , , 

ini \2T J 



duj , / LU \ 

coth — 

47ri \2TJ 



h ~ qi ki — Qi f de i9tanh(e/2T) 



(e - w)(e - uj')gR{k, e)gA{K e) 



47r de 

X [VAiq, w)VA(q', c^').9fl(k - Q, e - ^)gfl,(k - q', e - cj') - VA{q, w)Ffl(q', t^Os^k - q, e - io)gl(k - q', e - 
-l/fl(q, c^)FA(q', c^').9i(k - q, e - c^).g|^(k - q', e - + ^^(q, c^)Ffl(q', c^')ffi(k - q, e - c.).gi(k - q', e - a;')] 



(B5) 



r 



The above expression can be arranged in a rather com- 
pact form: 



de 9tanh(e/2r) 



d ^ I An 
k 



5fl(k,e)5^(k,e) (B6) 



I — coth 



dki ^-^ j Am 
q 



(^) (e-^) (VA(q,c.) 



X5K(k - q, e - tj) - V^Xq, c^)gA(k - q, e - t^; 



n 2 



Note that the expression in the squared brackets is reaL 
Without the factor e — cj, this expression is precisely 
the first order expansion with respect to the interaction 
of the renormalized velocity. For Eq. IB6I to be consis- 
tent with Eq. [63l the expression in the squared brackets 
must reduce into ed'Slai (k, e) /dki (where txi (k, e) is the 
self-energy term with a single interaction line shown in 

Fig.inra)). 

Below we shall examine the expression in the squared 
brackets following the calculation of the self-energy in 
the Fermi-liquid theory presented in Ref. Isil First, we 
rewrite this expression in the following way: 



A(k,6) 



d_ 

dk, 



n 



dqde' , f ^' \ i-,r , /\ 
- — —T-— tanh — ■ e Vr(ci, e — e ) 



^ n f dqduj / oj \ 

xXmgfl(k - O + y coth J 



X (e - uj)gR(k - q, e - uj)ImVR{q, uj) 



(B7) 



Next, we use the Lehman representation for the Green's 
function and the interaction propagator: 



A(k,e) 



_d_ 

dk, 



n 



dqde/ 

d+l 



i2n) 



duj 



e' tanh 



2T 



+ (e — uj) coth 



\2TJ 



ZmVflXq, uj)ImgR{k - q, e') 



LU + eJ 



iS 



(B8) 



One may approximate the imaginary part of the quasi- 
particle Green's function with a 5- function, Imgfl(k, e) ~ 
— 7r5(e — ^k)- For the sake of simplicity, we consider 
a three dimensional system with a parabolic spectrum, 
= fc^/2rn — /i. Then, the integration over the momen- 



tum q can be written as J dq = 2'Kq'^dq sinOdO: 

XmVR{q,uj) 



A(k,e) 



9 /■ de'duj 
dk, J J2nr 



q dq- 



LO - 



— e — iS 



e' tanh 
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2T 



(e — uj) coth { —=\ 
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fc2 - 2kqx + q^ 



m 



The integration over x yields non-zero result when (fc — 
qf/2m <e' < {k + qfl2m. 

Next, we transform from the integral diu to 

/o°° 



A(k,6) 



d 



(27r)3 dki k 



OO /'OO 

de'V I doj / qdq (BIO) 

^0 



^ f ^nT'VR{q,uj) 
\ LU + e'-e 
ImVR{q,u;) 



e'tanh( + (e - coth (^) 



e'tanhf — j - (e + c^) coth (^) 



where V J denotes the principle value. Finally, simple 
manipulations of the integrals yield: 

^(^'^)-(2^aA^Ty^^a '"/o 
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2uj coth I ) 
\2TJ 



tanh 
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(Bll) 



We may compare now between the expressions for the 
self-energy and A. In the self-energy the last two terms 
do not exist. The first term in the above equation for 
A is equal to the self-energy multiplied by the frequency 
e, and hence it precisely coincides with the renormalized 
heat current vertex (see Eq. [63]). Indeed, one can check 
that the integral over uj yields: 



Aa(k,e) 
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tanh 



27r)3 dh k 

e + e' 



2T 



j de' j qdqnVR{q,e') (B12) 

a7^(Tl(k, e) 
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2T 



dki 



This expression is exactly what we anticipate to get in 
order to satisfy the Wiedemann-Franz law. Therefore, in 
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the framework of the Fermi-liquid theory the contribu- 
tion from the other two terms in A must be zero. How- 
ever, only the last term is zero, while the second one is 
not: 



qdq (B13) 



ImVn{q,uj) 
uj + e' 
1 dm 



(27r)3 dh k 

e + e' 



e tanh 
de' 



e tanh 



e' tanh 



2T 



2T 

qdqnVR{q,e') 
e — e 



2T 





(' tanh 



2T 



One can see that this integral is not zero (the integrand 
is an even function of e'). Moreover, the integration over 
the frequency e' in Eqs . IB 1 2l and IB 13l accumulates over a 
range of frequencies that is restricted only by the typical 
scale of the interaction. We may conclude that the sec- 
ond term is comparable with the anticipated contribution 
of the first term (Eq. IB12P and, therefore, the difference 
between the renormalized heat current vertex and A can- 
not be neglected. Other second order contributions can 
not save the situation because in all of them the R-A 
section is connected to one of the bare vertices. Such 
contributions do not violate the Wiedemann-Franz law, 
see the discussion below Eq. lB3l where the first group was 
analyzed. 

In this appendix we calculated a particular contribu- 
tion to the thermal conductivity in the second order per- 
turbation theory with respect to the electron-electron 
interaction. We showed that already on the level of 
the Fcrmi-liquid theory, there is a disagreement between 
the quantum kinetic approach and the simplified Kubo 



formula. While the quantum kinetic approach repro- 
duces the known phenomenological behavior, the sim- 
plified Kubo formula fails this test. Therefore, the use of 
the simplified Kubo formula in the presence of electron- 
electron interactions is unjustified and may lead to erro- 
neous results. 



APPENDIX C: THE COULOMB DRAG 

In this Appendix we concentrate on the contribution 
of the Coulomb drag to the electric and heat currents. 
The Coulomb drag term is generated in the quantum 
kinetic approach when the external perturbation (either 
the electric field or the temperature gradient) enters the 
self-energy given in Eq. [55] through the propagator of the 
interaction. 

For the sake of simplicity, we will write the expression 
for the Coulomb drag term when the averaging over the 
disorder was already performed and therefore we can use 
the momentum space representation. The explicit deco- 
ration of this term with diffusons can be easily done after- 
wards. Therefore, we may Fourier transform the spatial 
coordinates in the expression for the self-energy given in 



Eq.[7i 



dqduj 



(Cl) 



(27r)'i+ 

-G<^>(k + q,e + c^)y>^<(q,w)] . 



Inserting the self-energy presented above into the third 
term in the expressions for the electric and heat currents, 
Eqs. [23land[57l one obtains: 



Je,h 



dkde 

1 f dkde dqdiu 



-vj'ik, 6)4(k, e)xe,h(e) [S<(q, e)(l - npie)) + I]>(q, e)npie)] (4(k, e) - ^(k, e)) + c, 



(C2) 



d+ 



jXeM^) Ki^^ e)fff,(k, e) - vfik, e)4(k, e)] [^^^^(k - q, e - - .g« (k - q, e - uj)] 



[4(k,e) -4(k,e)] [il-nF{e))nFie~u;)V<iq,iu)+nFie){nFie~u;) - l)V>i<i,u;)] 

'I J (2^)^+1 (2'^^)l+i ^^-^(') e)3fg(k, e) - vf{k, e)4(k, e)] [^(k + q, e + c^) - 5f,(k + q,e + u;)] 

[gf,(k, e) - 4(k, e)] [(1 - nF{e))nF{£ + uj)V> {q, lo) + nF(e) (71^(6 + ^) - 1)'^f (q, w)] 



Recall that the difference between the electric and heat 
currents is absorbed into the function Xe.h{^)^ where 
Xe(e) = — e while Xhi^) = e. In addition, F indicates 
the dependence on the external perturbation. Note that 
Eq. IC2I describes not only the electric and thermal con- 
ductivities, a and k, but also the two off-diagonal trans- 



port coefficients, a and a. One may connect the two 
contributions given above with the two possible locations 
of the velocity vertex inside the polarization operator. 
Using the identity for the products of the distribution 
functions given in Eq. I75bl the currents can be rewritten 
as: 
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Je,h 



dkde dqduj . ^ n r i?n ^ n 



(27r)'i+i (27r) 
[4(k, e) - 4(k, 6)] - xe,h(e - c^) [^;f (k - q, 



^(k,6)4(k,e)] [4(k-q,6 



4(k 



^)9eqO^ - q, e - - u/'(k - q, £ - w).g^g(k - q, e - u)] 



[5^9(k,e) -fff,(k,e)] [gf,(k - q, e - tj) - ^^^(k - q, e - w)] } [nj.(e) - n^Ce - w)] 



(C3) 



In the second term we have redefined the frequency 
e + uj ^ e. Here, U represents the two interaction fines 
connected to the polarization operator through which the 
external perturbation enters (see Fig. [9]): 

U = [np{-u;)V< + np{u;)V>]. (C4) 

Let us analyze the analytic structure of U. For 
this purpose, we need the expression for the F- 
dependent propagator of the interaction, VF(q, w) — 
-Veq{ci,uj)nF{q,uj)Veq{q,uj), see Eqs. [m and 111 The 
lesser and greater components of Vf are: 

(C5) 

As a consequence of the standard relations between the 
different components of the propagator at equilibrium, 
V< = np(^)[Kf - K^] and V> = -np{-u;)[Vj^ - V^], 
the non-vanishing contributions to Eq. IC4I occur when 
the two interaction propagators have an opposite analytic 
structure, V^V^. [If one identifies V^(q, uo) with the di- 
agram presented in Fig. fTST a) and V-^ (q, lo) with the di- 
agram in Fig. I13f b). it becomes clear that the symmetry 
with respect to q, u < — > — q — w is already embedded in 
lA. As shown in Ref. 29, this symmetry is responsible for 
vanishing of the Coulomb drag terms with two retarded 
or two advanced propagators.] 



Let us look at the lesser and greater compo- 
nents of the polarization operator, np'^(q, — 
/(ikde/(27r)''+iG<'>(k,e)G>'<(k-q,e- w). In the lin- 
ear response, one has to exploit the two possibilities to 
replace one of the Green's functions in IIf by the ap- 
propriate Gf presented in Eqs. [13] and [39l Examining 
Eq. IC3| one can recognize that the contribution for the 
electric and heat currents from the Coulomb drag can 
be presented in a very compact and symmetric way as a 
product of two polarization operators: 

[np{-uj)Tl<+np{Lo)I]>] [np{-uj)IV< + np{uj)Ii>] ; 



[np(-w)n^y -I- np(w)n^y] [np(-a))n^ -f np(a;)n^] . 

Next, we shall focus on one of the triangles in the 
Coulomb drag, [np(— w)np(q, -t- np(a;)np(q,Ci;)]. We 
separate the triangle into two groups of terms, Oi and 

<2- 



<ll = 



dkde dnpiu.) ^^^^^ [^B^i^^,^,)i^gRi^^,^,)fgA^^^^^^^^ 



-F(e - uj) [v«(k - q, e - c.)(.g« (k - q, 6 - u)f gf^{k, e) - v^(k - q, 6 - c.)(4(k - q, e - u:)f g^,{k, e)] } 

X [7ip(e - Lj) — np(e)]; (C7a) 



dkde dnp{uj) 



F(6){v^(k,e)(4(k,e)f4(k-q,6 



v-(k,e)(.9f,(k,6)f4(k-q,e-o.) 



(27r)'^+i duj 

[v^(k, e) + v^(k, e)] 5fg(k, e)4(k, e) [^(k - q, e - - g^^ik - q, e - c^)] } [np(e - w 



)-np(e)], (C7b) 



where F(e) is equal eE for the response to an electric 
current and eVT/T for the response to a temperature 
gradient. The expression for the second group, i.e., 
<l2, was simplified using the symmetry with respect to 



q, uj < — > — q, — cj. The analytic structure of this group is 
the same as obtained in the case of electric conductivity 
using the Matsubara techniquci^S We will show that the 
<! vanishes for the response to an electric field. Further- 
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(a) (b) 

FIG. 13: The two contributions to the Coulomb drag corre- 
sponding to W = [np{—uj)Vp{ci,uj)+np{uj)Vp{ci,Lo)]. (a) the 
diagrammatic interpretation for V^(q, cj) and (b) the diagram- 
matic interpretation for Vp{q,Lj) 



leading contribution from the Coulomb drag to the ther- 
mal conductivity has Oi on each side as illustrated in 
Fig. [9l This term is discussed at the end of Sec. |Vl 

To conclude, we showed that the analytic structure of 
the Coulomb drag terms is different for the responses to 
an electric field and a temperature gradient. In Fig. ll4l we 
present the corresponding drag diagrams. It is interesting 
to note that the combination <li V^V^ generated by the 
quantum kinetic approach cannot be obtained using the 
simplified Kubo formula described in Eq. We believe 
that the origin of this contribution is related to the part 
in Luttingcr's heat current operator— that depends on the 
interaction. 



more, we examine the fate of this term for the response 
to a temperature gradient. With this in mind, we re- 
formulate the products g^'^(k, e)v^''^(k, e)f;^''^(k, e) as 
a derivative of a single Green's function dg^^'^ik, e)/9k. 
After integrating by parts we obtain: 

Id f dkde t , . r a „ a •. 
'^^^2d^J ^^^'^ ^^^^'i^^ -^^('-^) 

^)9eq{^ - q, e - a;)] + F(e - ^) [.gf^(k, e) 

X5fq(k - q, e - tj) - 5^^^(k - q, e - t^)5e^,(k, e)] } 

X [nF(e-w) -nF(e)]. (C8) 

For the response to an electric field, when F(e) is inde- 
pendent of e, Eq. lCSl is identically zero. However, for the 
response to a temperature gradient, part of the RHS of 
Eg. IC8I survives, because F(e) and r(e — uo) are not the 
same. The remaining term is proportional to lo: 

Vr f dkde . , . , 
<i= -uj— J j^-:^[nF{e-Uj)-nF{e)] (C9) 

d 

0^ [fff,(k - q, e - ^^)gtqi^, e) - 9eqi^ - q, e - ^)geqiK e)] 

In the presence of disorder, this triangle can be decorated 
by two diffusons and, hence, its contribution to the cur- 
rents is significantly larger than the one from <l2- The 




(a) (b) 

FIG. 14: The analytical structure of the different contributions 
to the drag diagram for (a) the electric current and (b) the heat 
current as a response to either an electric field or a temperature 
gradient. 



APPENDIX D: EXAMPLES ILLUSTRATING 
THE ONSAGER RELATION 

In this appendix we present several examples in which 
we demonstrate the fulfilment of Onsager relation. We 
prove the Onsager relation for the last term in Eq. [8T] 
when the self-energy contains only one interaction line. 
This self-energy is described in Eq. [79] and illustrated in 
Fig. m Inserting the self-energy given in Eq. [79] into the 
last term in Eq. [52 we get: 

^V<{uj){l - npie)) + G<{e + Lo)V>{Lu){l - npie)) 
+G>{e - uj)V>{uj)nF{e) + G> {e + uj)V<{uj)nF{e)] 
(4W-.4W) + c-c. (Dl) 

Here either G = Gf is the field-dependent Green's func- 
tion and V — Veq is the equilibrium propagator or the 
other way around. 

We start with the case when the external perturbation 
enters through the Green's function of the quasiparticlcs. 
The F-dependent quasiparticle Green's function is given 
by Eq. [83] Using the identities for products of the dis- 
tribution functions given in Eq. [TS] we obtain that the 
contributions for the thermoelectric currents generated 
by Eq.[Dl]are 




FIG. 15: The diagrammatic representation of the contributions 
to the currents described in Eq. ID2I For simplicity we do not 
show here the dashed lines that indicate the scattering by the 
impurities. 
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[5f5(r5,r6,e)wj(r6,ri) (.g^^(ri, r^, e) - ^^^(ri, r^, e)) - (^^^(rs, r^, e) - gf^{r„re, e)) Vj{re,r,)g^g{r„r^, e)] 
[{9eqi^2,rs, e-uj) - 5fq(r2, r3, e - w)) Vj{r.„r^, )gf^(r^,r„e - uj) - 5f^(r2, Tj, e - uj)vj(r.„r^) 
(5^g(r4, r5, e - - 5fg(r4, r,, e - lu))] 

+ ^ 1 ^dr,...dr,^^^xe,h (e - (e) K (e - a.) - (e)] [y«(r, ,r,,u;)-V^ (r, , r, , u;)] 

[{9eq{^2,rs,e) -5fg(i-2,r3,e)) I'j(r3,r4)g4(i"4,r5,e) - .gf^Cr^, Tj, e)wj(r3, r^) (g4(r4,r5,e) - gfg{r^,r,, e))] 
[9eq{^5,r^, e - w)wj(r6, (.g^^ (r^ , , e - w) - gfg(ri, r^, e - t^)) - (^^^(rs, r^, e ~ uj) ~ g^^{r.^,re, e ~ w)) 
«j(i"6,ri)5fg(ri,r2,e- w)] • 

The two parts of the above equation originate from the two terms of the self-energy written in Eq. 1791 In the second 
part we have redefined the frequency e + w — > e. One can recognize that Eq. ID2I poses the structure of Eq. 1861 and 
as such satisfies the Onsager relation. Furthermore, we may reformulate the expressions for the electric current as a 
response to VT and the heat current generated by E in a compact way: 

^dr, {np{e~uj)~nF (e)] "MK? (r2 , , c.) - ^1 (r2 , rg , c.)] (D3a) 
X {G'^(r2,r5,e)G'^y(r2,r5,e -w) + G^j,(r2, rg, e)G^(r2, rg, e - w)] . 

VTj„ (E) = 1 ^dr, \np{t-uj)~ np (e)] (v^ ,r,,u:) - ^1 {^2 , rs , c.)] (D3b) 
X {GvT(i"2,r5,e)GE(i"2,r5,e - w) + GE(r2, rg, e)GvT(i"2, rs, e - cj)] . 

The diagrammatic representation of these contributions to the current is presented in Fig. 1151 It is worth nothing 
that one may identify A^qPeq appearing in the general proof of the Onsager relation (see Eq. [85]) with the effective 
interaction. Then, the contributions to the thermoelectric currents can be obtained from Eq. ID2I by substituting 
— with AeqPeq- Correspondingly, the compact form of the currents presented in Eq. ID3I can be reproduced. 
To explain what we have in mind, let us look at the self-energy given in Fig. [TUT b) and Eq. [Ml Following the same 
manipulations, the electric and heat currents generated by this self-energy can be written as: 

Eje(VT) J '^(2^)3^ dr^...dr^ [npie - uj - u') ~ npie)]^^ [np{uj) - np{-uj')] g^g(r,,r^,e - uj)gfq{r:,,r^, e - uj') 

{(Kg (ri, r3, w) - y4(ri, r3, c^))G< j,(r,, r3, e - u; - u'){Vj^{r,, r„ uj') - V^'^^ir^r,, ^'))G<{r„r„ e) 

-K(F^(r3, r„Lu') - V^ir,, v„uj'))G<j.ir„ v,,e^uj- uj'W{v,, r,,u;) - V^{r„ v,,Lu'))G<{r,,r„ e)] , (D4a) 

VTj,,(E) = -T^ J ^^^^dr,...dr, [npie ^ ^ - - npie)]'^ [np{uj) - np{-uj')] gf^{v,,T,,e - Lo)gf^{r„v,, e - J) 

{ ( (r 1 , r3 , c^) - f4 (r 1 , r3 , c^) ) G< (r, , r3 , e - ^ - c^O (Kg (r2 , , c^O - Kg (r2 , r4 , c^O ) G^T (^4 , r 1 , e) 

+ (F^(r3, r„ c.') - V^{v,,v„J))G<{v,, r„e-oj- u')(V''{r„ v,,w) ~ V^^(r„ r„ c.'))G< ^(r„ r3, e)] , (D4b) 

I 



The diagrammatic interpretation of these currents is pre- 
sented in Fig. [TH 

Let us return to the discussion of the self-energy pre- 
sented in Eq. [THl We shall examine the contributions 
to the thermoelectric currents in which the external per- 



turbation enters Sp through the propagator of the in- 
teractions. These contributions, corresponding to the 
Coulomb drag terms, were analyzed in Appendix [C] It 
is shown there that the electric current as a response to 
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f'i'fE-co-a)') 




F(e-o)-(n') 



FIG. 16: The diagrammatic representation of the contributions 
to the currents written in Egs. ID4al and lD4bl 



a temperature gradient can be written as: 

[np(-w)n^(r4, ri, w) + np(w)n^ (r4, v^,uj)\V^^{v^,Y2,uj) 
[np(-w)n^j,(r2, ra, w) + np(w)n^j,(r2, rj, w)] . (D5) 

Similarly, the heat current generated by an electric field 
can be formulated is 

[np(-cj)n^y (r4, ri, uj) + np(cj)n^y(r4, y^^uj)] (D6) 
V^{y,,v^,uj) [np(-w)n<(r2,r3,cj) +np(w)n>(r2,r3,w)] . 



One may see that Eqs. lD5l and lD6l are connected through 
the Onsager relation under the condition of microscopic 
reversibility. However, it is not obvious how to recognize 
in the above expression the structure of Eq. [M] used in 
the general proof of the Onsager relation. To resolve this 
issue, one has to go to Eq. IC2I in Appendix [Cl Then, 
instead of analyzing each polarization operator IIf sep- 
arately, one should concentrate on combinations of the 
kind [(1 — np{e))np{e — aj)np(q, oj) + np{€){np{e — to) — 
l)np(q, w)]. Following Sec. IVIi we shift the arguments 



of the Green's functions inside SF(e) in such a way that 
they all contain e. The polarization operator will be writ- 
ten as: 

n<'>M = i|^[G<->(.-l])G>-<(e-17-c.) 

+G<^> (e -I- 17)G>'< [e + VL^uj)]. (D7) 

Now, the Coulomb drag term contains four distribution 
functions. To simplify the expression, we group these 
functions into two pairs, and use the identities in Eq. l75bl 
for each of the pairs. Here, we pair the distributions func- 
tions in such a way that the difference between the two 
frequencies in each pair is O. Namely, one pair contains 
np{e) and np(e — 51), while the other includes np(e — w) 
and np(e— w— rj). Following the steps described here, one 
can obtain an expression for the Coulomb drag term that 
has the same structure as Eq. 1861 Consequently, one may 
present the Coulomb drag contributions to the thermo- 
electric currents in a compact form similar to Eqs. ID3al 
and [D3bl 

APPENDIX E: THE QUANTUM KINETIC 
APPROACH IN THE PRESENCE OF 
SUPERCONDUCTING FLUCTUATIONS 

In this work, we derived the currents as a response 
to an electric field and a temperature gradient for elec- 
trons interacting through the density channel, e.g., the 
Coulomb interaction. Here we show that following the 
same scheme one can find the expression for the currents 
in the presence of superconducting fluctuations. 

The difference between the fluctuations in the two 
channels reveals itself most clearly when the electric cur- 
rent or the response to an electric field are studied. The 
action for electrons interacting through the Cooper chan- 
nel in the presence of an electric field is 



A(r, (r, (r, t) + /I.C 



(V^t(r,t))(VV,^(r,^)) 
2m 



^ Vl(r, t) [er • E + F„„p - ^] ^.(r, t) (El) 



A(r,t)| 



It follows from the action that the charge of the 
electrons, — e|'0(r)p, is not conserved unless the 
charge carried by the interaction field, 2e|A(r)p, 
is also included. Correspondingly, the continuity 
equation for the current of the quasiparticles ac- 
quires a source term — e c^t IV's(r, i)P + Vj^P(r,t) = 
-2ie[At(r,t)V^;(r,t)VT(r,t) - A(r, t)V| (r, t)^! (r, i)]. 
Hence, the expression for the electric current is the sum 



of the currents of the quasiparticles and interacting field: 



"2 



dr'dt' 



2v(r,t;r',t')G(r',t';r,t) 



(E2) 



+ h.c. 



Here, L^{r,t;r',t') is the gauge invariant propagator of 
the interaction in the Cooper channel, and we use the 
notation V(r, i; r', t') = —i{r ~ r')tlg^{r,t;r' ,t') to em- 
phasize that this term is analogous to the renormalized 



29 



velocity given in Eq. [Ml [Since the propagator L does 
not have a dimension of an inverse energy, V does not 
have the dimension of velocity] The factor 2 in the con- 
tribution to the current from the Green's function of the 
quasiparticles is due to the sum over the two spin direc- 
tions, while the factor of 2 in the contribution of L is 
because the superconducting fluctuations carry a charge 
of 2e. 

Since A is a charged field, the kinetic equation for 
its propagator L(r, t; r', t') resembles the kinetic equation 
of the quasiparticle Green's function, see Eq. [TUl rather 
than Eq. [T5]for the neutral interaction field: 

— X^^L{p, uj; imp) — 6{p) — / dvifl (p — ri, cj; imp) 



L (ri, (jj; imp) . 



(E3) 



Using Eqs. IE2I and IE3[ we have reproduced the known 
results for the paraconductivit y'^^i'^^'^^ and for the cor- 
rections to the magneto-resistance in disordered super- 
conducting films. '^^ 

Let us turn to the response to the temperature gradi- 
ent, which is analyzed with the help of the gravitational 
field. Unlike the response to an electric field, the quan- 



tum kinetic equation and L in the presence of a gravita- 
tional field has the same form as Eq. [4T]for V: 

- A-iL(r, t- r', t') = 5{r - r')5{t - t') (E4) 



dv^dtjl (r, t; Y^,t^)L{Yi,t^;Y\ t'). 



Here, L_ is transformed according to Eq. [29l Note that 
just like for y, the dependence of the propagator L on 
the gravitational field is only through its self-energy 11^^. 
In view of the similarity of the kinetic equations for V_ 

and L, it is natural that the heat current (both as a re- 
sponse to E and VT) does not acquire any additional 
terms in the presence of fluctuations in the Cooper chan- 
nel and Eqs. [HT] and [55] are still valid. Finally, the ex- 
pression for the electric current as a response to a tem- 
perature gradient includes two contributions; one is the 
electric current of the quasiparticles (described by the 
VT-dependent Green's function) and the other is from 
the electric current carried by the interaction field (de- 
scribed by the VT-dependent propagator L). 

We applied this scheme for the calculation of the 
Nernst effect in disordered films above the superconduct- 
ing transitionJ-5. 
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